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ABSTRACT 


This document contains the Spectral Analysis Program (SAP) flow 
diagrams and listings. The Plot Generation Program (PLTGEN) flow diagrams 
and listings are also included. 
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1. INTRODUCTION 

The purpose of this nianual is to present the listings and flow 
diagrams for the Spectral Analysis Program (SAP) and the Plot Generation 
Program (PLTGEN). SAP computes the power spectrum of an angle modulated 
signal and is useful in the analysis, design, and testing of angle modulation 
communication systems. PLTGEN accepts a SAP output tape and generates the 
appropriate commands required to drive the EAI Dataplotter. The PLTGEN 
output tape is mounted on the Dataplotter to obtain a power spectrum plot. 

The detailed mathematical formulation of SAP is presented in Reference 
1 (Final Project Report). A user's guide, providing the necessary information 
to 1) understand the general SAP computational approach and software struc- 
ture, 2) set-up the input parameters in the appropriate tape or card format, 
and 3) execute the program, is presented in Reference 2 (The Spectral 
Analysis Program, Volume I - User's Guide). User information for PLTGEN 
is also included in Reference 2. 

SAP and PLTGEN were written using the standard TRW self-documentation 
technique. With this technique, each subroutine contains: 

• Programmei -nd date 

• Purpose of subroutine 

t Description of each input/output parameter on cards, tape, 
or through COMMON 

• Remarks and restrictions 

• Additional subroutines required 

• Numerous functional descriptions of sections of programming 
logic. 
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2. SPECTRAL ANALYSIS PROGRAM 


V 


c 




This section presents 1) a brief description, 2) flow diagrams, and 
3) listing for each module or subroutine used in the Spectral Analysis 
Program. The presentations are in alphabetical order. 


2.1 SAP MODULE/SUBROUTINE DESCRIPTIONS 

The following list describes all of the modules/subroutines for' the 
SAP program. 

BITREV Bit reversal routine 

CL0CK References computer clock 

EFFT Computes extended fast Fourier transform 

FFT Computes fast Fourier transform 

. FILTER Performs filtering of the modulating or modulated 

signal 

ISAR Generates a signal tape compatible with SAP from 

a user supplied input tape 

MAGTAP Tape operations routine 

M0D Performs the exponentiation operation used to 

generate the modulated signal representation 
from the modulating signal 

PL0T Generates the SAP plot tape 

SAP Spectral Analysis Program driver program 

TIMBl Computes timing differences 

TRFN Calculates the real and imaginary parts of the 

desired transfer function 

TSGEN Generates built-in modulating test signals 

TTRANS Large matrix transpose routine 

TWRITE Prints signal at specified stages of processing. 
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2.2 SAP FLOW DIAGRAMS 

A complete flow diagram for each of the module/subroutines listed 
in Section 2.1 is given in the following diagrams. 










DUMMY TIMING ROUTINE* (CLOCK) 



* The sole purpose of this dummy routine is to allocate 
enough core storage for a CS-1 routine named CLOCK 
that will be loaded later from paper tape. 



Figure 2-2. Dumny Timing Routine* (Clock) 
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Extended FFT Subroutine 


ENTER WITH 

UNT1 

UNT2 

NR 

NC 

\ 

unit number of input tape 
unit number of output tape 

number of rows in input array (must be a power of 2) 
number of columns in input array (must be a power 
of 2) 

MODE } = 

1 compute direct FFT 
f compute inverse FFT 

AR 


AI _ 

WR 

working arrays of dimension max (NR,NC) 

MI j 



Rewind TaO' 
UNT1 & UNT 

a Units 

i 


Determine 1 
NR = 2 N+ 

'i such that 
1 


SET l 

1 

-2tt 

IK 

if MODE « 1 

C0N*=j 

% 

2* 

NftxWC 

If MODE i 1 


Figure 2-3, Extended FFT Subroutine 







m 



gure 2-3. Extended FFT Subroutine (Cont’d) 






c 



% 



Determining N such that 
NC * 2* N+] ) 


FFT each column of array on UNT1 


FOR J = 1 , NR 

1) Read next two records on taoe unit UNTI and store 

in AR(1 ) and A1(I), I = 1, , NC 

2) Set NCW = J 

3) CALL FFT (N, MODE, NCW) 

4) Wri te AR and AI as next two records on tape UNT2 


Rewind tape units 
UNTI , UNT2 


EXIT 

Resulting complex array (NC rows and NR columns) 
Is stored by columns on tape unit UNT2, N. B. 
independent variable increases along a row. 




Figure 2-3. Extended FFT Subroutine (Cont'd) 
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Fast Fourier Transform Routine (FFT) 


C 


( 



<b 


Figure 2-4. Fast Fourier Transform Routine (FFT) 
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FILTER MODULE 


ENTER WITH 


I TAPI - Input tape number 
ITAP2 = Output tape number 

NC = Number of columns in ARRAYS 
NR = Number of rows in ARRAYS 
DFREQ = Frequency increment 
FC = Carrier frequency 

! 1 Butterworth 
2 Chebychev 
3 Input Data Filter 
4 Integration Filter 
5 Product of 4 & 1 
6 Product of 4 & 2 
7 Product of 4 & 3 

TMAG 

TPH Data as required by 
F transfer function routine 
NF 

NORDER = ORDER 

FCEN = Center of frequency 
FCUT = Cutoff frequency 
RIP = Ripple 


Rewind tape 
Units ITAP1 
& ITAP2 

SET I FLAG & I OPT BY 

IFIL 1 2 3 4 5 6 7 

I OPT 12 3 12 3 

I FLAG 1 1 1 2 3 3 3 



Figure 2-5. Filter Module 
2-12 














Read next two records from 
Tape ITAP1 , each record has 
NR words. Store as 

1st Record AR(I) 1=1, NR 

2nd Record AI(I) 1=1, NR 


I FLAG 
= 2 



Filter AR and AI for 
frequencies above FC 

1. For 

1 = 1, NR/2 

a) 

Set 


FREQ = FLOAT ((0-1) + (I-1)*NC)*DFREQ + FC 

b) 

Call TRFN with FREQ, I0PT, PAR, TMAG, 
F, NF, and return with TR and TI 

c) 

Set 

ARP = AR (I) 
AIP = AI (I) 

d) 

AR (I) = ARP*TR - AIP*TI 
AI (I) = ARP*TI + AIP*TR 


Figure 2-5. Filter Module (Cont'd) 










Filter AR and AI for frequencies 
below FC (Modulated Signal) 


1. For I = (NR/2) + 1, NR 

a) Set IP = NR + (NR/2) +1-1 

b) Set 

FREQ = FC - FL0AT(NR*NC-(J-1 ) 

- (IP-1 )*NC)*DFREQ 

c) Call TRFN with FREQ, I OPT,. PAR, TMAG, 
F, NF and return with TR and TI 

d) Set 


ARP = AR(IP) 
AIP = AI ( IP) 


e) Set 


AR(IP) = ARP*TR - AIP*TI 
AI(IP) = ARP*TI + AIP*TR 


m 
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Integration filter - modulating 
signal only 

1 . For 1 = 1,’ . . 0 

a) FREQ = F! { (0-1 ) 

+ u-i)*nc)*dfreq 

b) TI = -1 ./FREQ 

c) I F( FREQ . LT . 1 . 0E- 30) TI = 0 

d) ARP = AR(I) 

AIP = AI(I) 

e) AR(I ) = -AIP*TI 
AI(I) = ARP*TI 

2. For I = (NR/2) + 1 , NR 

a) FREQ = FLOAT (NC*NR - (J-l) 

- (I-1)*NC)*DFREQ 

b) TI = +1 ./FREQ 

c) I F( FREQ. LT. 1 .0E-30) TI = 0 

d) ARP = AR(I) 

AIP = AI (I) 

e) AR(I) = -AIP*TI 
AI ( I ) = ARP*TI 


m 









INPUT SIGNAL MODULE 


Enter with I 

IUNT1 = Input tape urit 

1UNT2 = Output tape '-nit 

NRI = Number of rows in input array 

NCI = Number of columns in input array 

DTI = Sampling interval for input data 

NR2 = Number of rows in output arrays 

NC2 = Number of columns in output 

arrays 

KSR = Defines output samplinq interval 
as DTI /KSR 


Rewind tape units 
IUNT1 & IUNT2 


Read record from unit IUNT1 
and store in S(I) I * 1 , NRI 


Figure 2>6. Input Signal Module 








Initialize 

AI(I) ■ 

0. 1 = 1, NR2 

PI = 

s(i) 

P2 * 

S (2 ) 

11 = 

2 

12 = 

1 

J1 * 

1 

02 = 

1 

TERM = 

( P2 - P 1 ) / FLOAT (KSR) 


Interpolate for Output Point 


AR( 1 2 ) = ( Pi +TERM* FLOAT (K-l )) 
12 = 12 + 1 
K = K + 1 


12 > NR2 



Write two records on unit 
IUNT2, each has NR2 words 

1st AR( I) 1=1, NR2 
2nd AI(I) 1*1, NR2 


Figure 2-6. Input Signal Module (Cont'd) 
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Figure 2-6. Input Signal Module (Cont'd) 






Dummy Tape Control Routine* (MAGTAP) 


ENTER WITH 


Parameter 

■j - Dummy Call 

i 



A(I) 

= l.o| 

J 



r > 

RETURN 


The sole puroose of this dummy routine is to allocate enough 
core storage for a CS-1 routine named MAGTAP that will be 
loaded later from paper tape. 


Figure 2-7. Dunmy Tape Control Routine* (MAGTAP) 



Modulation Module (MOD) 


ENTER WITH I 

IUNT1 = input tape unit 
IUNT2 = Output tape unit 
NR = Number of rows in 
input arrays 

NC = Number of columns in 
input arrays 

BETA = Scaling parameter 


Rewind IUNT1 
and IUNT2 


Compute exponential for each term 


FOR 0 = 1 , NC 

1) Read two records from IUNT1 and store in 
AR(I) , AI(I), I = 1, NR 

2) For 1=1, NR 

AR(I) = COS (BETA*AR(I )) 

AI (I ) = SIN (BETA*AI(I» 

3) Write two records on IUNT2 

AR(I) , 1=1, NR 
AI ( I ) , 1 = 1, NR 


Rewind IUNT1 
and IUNT2 


Figure 2-8. Modulation Module (MOD) 
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PLOT MODULE 


C 


Enter with 


1UNT1 

= 

InDut data tape unit 

IUNT2 

= 

Outnut tape unit for plotting 



equinment 

NR 

= 

Number of rows in inout ARRAYS I 

NC 

= 

Number of columns in input ARRAYS 


( 

> 

0 Frequency (in hertz) at which 


1 


Dlot disDla.y begins 

FSTRT 



- ZOOM type disnlay 



< 

0 Carrier frequency in center 


V 


tyoe of display 

FC 

= 

Carrier frequency (in hertz) 

NAVE 

= 

Number of Doints to be averaged 



for each outDut Doint 

DELF 

= 

Increment in frequency variable 



between input Doints 

AMAX 

= 

Maximum nower component of input 



signal 



If .LE. zero module computes and 



outputs AMAX 

IPMODE 

= 

0 

Print and Dlot 



1, 

Plot only 



2, 

Print only 

IPLPOS 

= 

0 , 

Print on unoer half of clot 




caper 



1, 

Print on lower half of plot 




paper 


Rewind tape units 
IUNT1 and IUNT2 



( 


Figure 2-9. Plot Module 
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COMPUTE MAXIMUM POWER COMPONENT 


For each J = 1 , NC 

1) Read two records from unit 
IUNT1 , each has NR words 
and store in 

AR(I) 1=1 NR 

AI(I) 1*1,. ..,NR 

2) Determine maximum of AMAX 
and (AR(I)**2 + AI(I)**2) 
1=1,..., NR and store in AMAX 


Rewind tane unit 
IUNT1 


NPS = 480 


! Determine Plot Tvne 

Set 




IFC * 0 

IF 

FSTRT < 

0 

IFC = 1 

IF 

FSTRT ^ 

FC 

IFC « -1 

IF 

FSTRT < 

FC 


i rnoo “ i 

NPS - NPS/2 


IFC * 0 


Figure 2-9. Plot Module (Cont'd) 










XI 


Check bounds and adjust NAVE if 
necessary 


1) Set FMAX * DELF*FL0AT(NR*(NC/2)-l) 

FLAST * DELF*FL0AT(NAVE*(NPS-1)) 

2) If FLAST FMAX then set 

NAVE = ((NR*NC/2)-l)/(NPS-l) 

and write out "NAVE has been adjusted 
to (value of NAVE)" 


Check bounds and adjust NPS if necessary 


1) FLAST = FSTRT + DELF * FL0AT(NAVE*(NPS-1 )) 

2) If IFC = 1 Set 

FMAX = FC + DELF * FLOAT ( (NR*NC/2) -1 ) 

If IFC = -1 Set 
FMAX = FC 

3) . If FLAST FMAX then set 

NPS = ( I NT ( (FMAX- FSTRT)/ DEL F)+l) /NAVE 

and nrint out message "number of ocints 
dotted has been adjusted to (value of NPS)" 


Initialize 


EPS * l.E-7 
SUM * 0. 
NSUM » 0 


PMIN « -70. 

SAVE - FLOAT (NAVE) 


Figure 2-9. Plot Module (Cont'd) 
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I PASS = I PASS + 1 




Print out P(K) and 
F(K) K = 1, NPS 


PMODE 
= 2 


Reinitialize for lower half of FC 
centered plot 


SUM = 0. 

NSUM = 0 
K = 1 

NSTART = NR * NC + 1 - NPS * NAVE 
N = NSTART/NR 
KT = NSTART - N * NR 
NSTRT = FC - DELF * FLOAT (NPS*NAVE 


Rewind tane 
unit IUNT1 


I General Tane for Plotting Equipment, 1 

1 Call 

PLTGEN 

a) 

Tane unit is IUNT2 

b) 

P(K) , F(K) , K = 1 NPS, are the 

dependent and independent variables. 


respectively 

c) 

P(K) is to be plotted as a vertical 
line for each K 


Figure 2-9. Plot Module (Cont'd) 















SCALE SIGNAL 


1) Rewind Unit IDT 

2) Set SMAX = 0. 

3) For J=1 , NCIDT 

a) Read two records from IDT 
Store in 

AR(I ) 1 = 1, NRIDT 

AI(I) 1 = 1 , NRIDT 

b) Set 

SMAX = MAX [SMAX, |AR(I)| 1=1, NRIDT] 

4) Rewind IDT 

5) For J=l, NCIDT 

a) Read two records and 
Store in 

AR( I ) 1 = 1, NRIDT 

Aid) 1=1, NRIDT 


b) Set 

AR(I) = (l./SMAX)*AR(I ) 1=1, NRIDT 

c) Write two records, AR(I) and Aid) 

1=1, NRIDT, on Unit 1ST 

6) Rewind IDT and 1ST 


Figure 2-10. Main Program (SAP) (Cont'd) 







*Ivc and Ivr are special notations, for convenience only , 
defining the relationship between the independent variable 
(time or frequency) and the rows and columns of the arrays 
at the point indicated as follows: Ivc - denotes that the 

Independent variable increments along a column (or within 
a record), thus, the next noint is generally in the same 
record while Ivr - denotes that the independent variable 
increments along a row (or from record to record), thus, 
the next point Is generally in the next record. 


Figure 2-10. Main Program (SAP) (Cont'd) 
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Yes __ 

(Filter modulatlnq 
Signal) 

(Integrate signal) 



Figure 2-10. Main Program (SAP) (Cont'd) 



































Ti mi nq Difference ComDUtation Routine (TIMER) 


ENTER WITH 

all to computer clock 
call to computer clock 

I = First c 
J = Second 

— 
j 


Compute Difference in Seconds 


X = J = I 
T = X/ 1 024. 

M = T/60. 

Y = M 

S = T - Y*60. 


RETURN 


Figure 2-11. Timing Difference Computation Routine (TIMER) 











A 



Figure 2-12. Transfer Function Computation Routine (Cont'd) 
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TEST SIGNAL GENERATION MODULE 


r 

V. 


Enter with 


IUNT = Tane unit for output signal 
NR e Number of rows in outout arrays 
NC = Number of columns in output arrays 
DT = Delta time between output points 
1 [Sinusoidal test signal 

I0PT = 2 (Periodic four level grav signal test signal 
3 (Periodic square wave test signal 

A(I) ) Amplitude, frequency, and 
F( I ) ( phase for IOPT * 1 
e(l) i option, each has 
M J M ^ 25 components 

TP (1-4) | For ICPT = 2, defines time points 
A(l-4) ) of breaks and change in amplitude 
) at breaks 

TP(5) Defines period of four level 
test signal 

TP(1), I For IOPT = 3, defines period and 
A(l) j amplitude of square wave 


Rewind tape 
unit IUNT 



c 


Figure 2-13. Test Signal Generation Module 
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I OPT * 1 


For I = 1 , NR 
M 

a) AR(I) * A(L) * 

L=1 

SIN(F(L)*T+e(L)) 

b) T = (I-1)*DT + TSUM 


Figure 2-13. Test Signal Generation Module (Cont'd) 














Figure 2-13. Test Signal Generation Nodule (Cont'd) 
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Transpose Subroutine 


ENTER WITH 

UNTIN 

UNTOUT 

NR 

| NC 

= unit number of input taoe 
= unit number of output taoe 
= number of rows in input array on UNTIN 

= number of columns in input array on UNTIN 

1 


Transpose complex array on tape unit UNT1 and store 
result on tape unit UNT2. Arrays are stored by columns 
alternating columns being the real and imaginary com- 
ponents. BUF working area is to be used to minimize 
the number of rewinds required of input tape. 


i 


EXIT 1 

The transposed complex array on tape 
UNT2 has NC rows and NR columns 


Figure 2-14. Transpose Subroutine 
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Taoe Signal Read and Print Routine (TWRITE) 


ENTER WITH 


IU = Data tape unit 

NC = Number of record pairs 
on tape 

NR = Number of words/record 


Rewind IU 


Read data tape and print data 


FOR J = 1, NC 

1) Read next two records on taoe IU and store in 
AR(I ) and AI(I), I * 1, .... NR 

2) Print J 

3) Print real part of siqnal ARd). 1*1, NR 

4) Print imaginary part of siqnal Aid), 1=1 . NR 


Rewind IU 


RETURN 


Figure 2-15. Tape Signal Read and Print Routine (TWRITE) 










2.3 SAP LISTINGS 

This section presents a complete listing for the ShP program. 


SAP PkObKAM LXSTIuG 


fiAXUpKOcKM. jkP 

Iw ♦*»♦♦■ ♦♦♦*♦*♦♦****♦****♦*♦*******♦*♦ 

CD . 

Lu KkooK*;-,. .cR ’\.jO tATu 

CU KXC.if.KU u » IftOl-iAS 

lu ; Tt<»t SlbTt.'.S 

UU / UUNt 1 ll<L 

uu 

uL POuFoji. 

Co Ti.xs xb 1 iic nAiri KOUTInt* IT REaDS ALL USER INPUTS AND 

Cl. c^ntr^jub l r it. PuOw through the subroutine 

OL. 

tL otiC'KXr' l XO.< Or PA,\k,*,E TE«S 

Co 

Co X.»Pol 

Cl. 

CO Ckkl^ 


tu 

iwf 

- 

uATA 

(APE UNIT 


C u 

leP 

- 

Pr.GoK«M ENTRY POiNT 


CU 

IP-1 

- 

• c.0* 

1 

FM MGuULmTION 


cu 


• 

• tut* 

2 

phase modulation 


CU 

IFOPl 

- 

• CO • 

0 

NO FILTERING 


CU 


- 

• L.U • 

1 

filter modulating of test signal 

ONLY 

Cl/ 


• 

• LU • 

2 

FILTER MODULATtD SIGNAL ONLY 


CU 

• 

- 

• LU • 

3 

FILTER TEST AND MODULATED SIGNAL 


C U 

xr CYPi 

- 

• Lu • 

1 

dUTTLRWOKTn FILTER 


tLf . 


- 

• Lvj • 

2 

CHEdYCHEV FILTER 


CJ 


- 

• L»tf • 

3 

INPUT FXLTlK 


cu 

xPTYPe 

- 

• CU* 

1 

GUI TLRwORTh filter 


cu 


- 

• tu# 

2 

CHtbYCHEV FILTER 


cu 


- 

• tv* • 

3 

INPUT FIlTuK 


U) 

xPLPoj 

- 

• cu* 

u 

PLOT «ILL BE ON UPPER HALF OF PLOT 

PAPER 

Cu 



• Lu • 

1 

PLuT WILL aE ON LOWER HALF OF PLOT 

PAPF 

Cu 

XP.'iOut 

- 

• tu* 

u 

PRINT POwCK SPECTRUM AND GENERATE DATA 

cU 



• cu • 

1 

ONLY GENERATE DATA TAPE 


cu 



• Lu • 

2 

ONLY PRINT POWER SPECTRUM 


cu 

xawp 1 

- 

• Lu • 

0 

ustR input signal on tape idt 


Cu 


- 

• LU • 

1 

SINUSOIDAL TEST SIGNAL 


cu 


- 

• Lu* 

2 

PEkIOuIC Four level GRAY SIGNAL TEST 

cu 


- 

• Lu* 

3 

periodic square wave test signal 
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O/ 

cu 

Lb 

Lb 

Lb 

Lb 

Ll^ 

u 

Lb 

bL 

Lb 

CL 

tu 

Ci) 

tb 

Cb 

Cb 

Lb 

bb 

Cb 

bb 

Lb 

Cb 

cu 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 

Cb 



SAP PROGRAM LIST InG 
IjI “ SPhRE TAPE U.-Ill 

K-»n “ iNc.3 SrtinPLluG RATE FOR RESAMPLED SIGNAL 
NoUF - bi*L OF bUFFER FuR TAPE TRANSPOSE ROUTINE* TTRANS 
l.blb'l - imJinuEk Op COLUMNS I j 4 bATA MATRIX 
url - wtFlhcS nUMbt-K OF PblNIS FOR INFUT FILTER* 
iFlYPl -• 3 

NPZ - ubFINES .MUiibLK OF POINTS FOR INPUT FILTER* 
lF 1 YP*: - 3 

" • hjnoEK OF KOaS Ik UaTA MATRIX 
I.Ai - mUV.jER of RECORDS On INPUT LATA TAPE 

NAVE - Nb,.uEK OP POIuTS TO SE AVERAGED IN OUTPUT 

ItPrtT - i Hi S mRRAY CONTROLS THE PRINTING OF THE SIGNAL 

o.< ‘fHL LATA TAPE. THt SIGNAL ON THE TAPE WILL dF 

pKAuiob br a call to twrite if the parameter is 

luomL TO 1* 

AIHkIca) nt-FtRS TO print after tsgen or ISAR 

ICRKIlt.) nc.Fc.KS TO PR ini AFTER FiRST CALL TO EFFT 

UPKTLj) INFERS TO PRInT AFTER FIRST CALL TO FILTER 

I<-PRU*t) kLFLRS TO PRINT AFTER SECOND CALL TO EFFT 

UPRllS) ntFERS 10 PKI..T AFTER TuIRO CALL TO EFFT 

UFKHd) litFEKS TO PRINT AFTER SECOND CALL To FILTER 
Ai.»mX - MAXlHbM POWER COr-.PONENT OF INPUT SIGNAL 
bET* - iiOOULmTIOn iNbtX. MULTIPLIES ENTIRE TEST SIGNAL 
bEbT - Hint SEPAKATIUNC IN SECONbS) bETWEEN POINTS IN 
ItST SIGNAL OR BETWEEN USER INPUT SIGNAL PRIOR 
lO RESAMPLING 

1 1 1 Ti-w - UP rc 30 CHARACTERS FOR RUN TITLE ANU/OR DATE 
FC - CARRIER frequency 

Hi - rNEuUcNCYL IN HER! 2) AT WHICH MAGNITUDE AND PHASE 
IS SPtCIFIEU FOR modulated SIGNAL 
F* - rKEGUENCYCIu HERTZ) AT WHICH MAGNITUDE AND PHASE 
IS SPtCIFIEU FOR MODULATING SIGNAL 

fceni - bturtn frequency for butterworth or chebychev 
fcen* - center frequency for bi .terworth or chebychev 
R lPi - RlPPLc. for butterworth filter 
r.iAGi - magnitude of input filter transfer function 

Tin AGE - WA3NI1U0E OF InPuT FiLTER TRANSFER FUNCTION 
NLKUi - uRLER OF THE bUT TERWORTH OR CHEBYCHEV FILTER 
NbRbZ - ORuEK OF THE bUTTERwORTH OR CHEBYCHEV FILTER 
FCJU - CUfoFF FREQUENCY FUR BUTTERWORTH OR CHEBYCHEV 
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CO 

PCUl'R 

- 

Curol-F FREQUENCY FOR BuTTER/.ORTH OR 

CHEBYCHEV 

cu 

KiP^ 

- 

RIPPLc. FoK ouTTEkWORTH FILTER 



CL* 

li'Hi 

- 

CORRESPONDING PH»\SEL IN KALIANS) 

- 


Cu 

Ipml 

- 

CORRESPONDING PHmSECIN PAOIArS) 



CU 

w 

- 

i»Ui«',t»El\ UF SlWUSOxuSL LL.ES) 



CL* 

AC i j 

- 

Mt-i^LlTUOE OF I SINUSOlu 



CC 

Cl I ) 

- 

FKEOUENCYLIN HlRIZ) OF I SINUSOID 



Cu 

T nc. T A 

- 

PMhSl ANGlELIn UtOKtES) OF I SINUSOID 


CU 

A 

- 

• luPT .EQ.id CHmNGc. IN AMPLITUDE OF BREAKS 


Cu 

M 

- 

•Iupt.eq.3 amplitude of square wave 



Cb 

IP 

- 

.XOF’T.EO.3 PERIOD uF SQUARE WAVE 



CU 

TmPu 





cu 






cu 

uulpul 





Cu 

Cu.l*.Oi4 





CU 

PranT 





cu 

LuCAT 

- 

LHSI ENTRY POINT PASSED 



cu 

IuT 

«* 

t »Ta tape unit 



Cu 

1ST 

- 

*PARE TAPE UNIT 



cu 

NKluT 

- 

NUMBER OF RECORD PAIRS CURRENTLY ON 

DATA 

tape 

;o 

NCXuT 

- 

number of worqs/recoro currently on 

41T1 

T175 

cu 

WIN 

m 

TOGETHER WIN AND SECS GIVE THE ELAPSED TIME OF 

cu 

St.CS 


SEGMENT OF THt PROGRAM. 



cu 

UPt. 





cu 






cu 

kEnAKKS ANu RE j 

IK I CT I ON* 



cu • 

GUahuul ii»ES REQUIRED 



cu 

BiTRtV 





cu 

ClOCis 





cu 

EFFT 





cu 

FFT 





cu 

FiuTEK 





CD 

lSAR 





cu 

MhGTAP 





CU 

KOu 





Cu 

PLOT 





CU 

TirtER 





cu 

TkFN 





cu 

TSGEN 





cu 

TTKANS 
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Co T ArUTc 

cu 

CU* *************#********. *****4****** 

oiMLnSiu.M ZlC2U4tt) rZ2L2048) 

uii'it.oiuN aklZdo) tAitdSS) r wKL25b) »wK2S6) 

uiMtuSiGW FlCbu) »F«:LbO) »PARUbO) t PAR2L 50 ) t TMAGlt. 50 ) #TMAG2C50) 
i » 1 PhiL ,»o ) » i PH2L bO ) »AL<ib) *FL2b) »T h£TAC 2S) »TPC6) » ITITLEC 30) 
ul.'itu IZPKfLb) 
uI.mLi^sagn imaglcS) »rpntzb) 

COi'iMuW/bTOrto£/Ar<» Al r 21*Z2 

COMmoN/UTu/ITITLE 

c^Oi vAuci^Cc-L WKC 1) f /iL i) ) 

LviOi t/A»_c_ijCLL WiC 1 ) » Z2C 1 ) ) 

CALL CuGCKL IT IMLI ) 

11-1 
iZ = c 
l4 = 4 
17 = 7 

CALL KhoI API i.4 » 14 f AR t II ) 

CALL M AO T Apt. 1 4 1 1 4 » AN f 1 2 ) 
hcaug»iuuz) UTirLtci) #i=i» 3 o) • 

1UU2 FORMAT CdUAl) 

KlaOl br 1001) IoT * IlP #IFM rlFOPT » IFTYP1 » IFTYP2» IPLPOS* 

A iPMOudrlSOPT rIST »KSR #mC1UT »NC1 »NF1 »NF2 » 

2 ,ma 1E»T fUKl »NAVL #CIZPRTC1) »I=1»6) 

1UU1 FORMAT! ibid) 

. htAui. G» 1003) amaX»u£TA»D£LT»FC»FSTrT 
cRi TtLo* 199d) 

199d FOKMATLlril) 

wRITiiCo* 1999) 

1999 FOKMATL 1H t lOdd*************************************************** 

4 ******* ************* ********************** ******** **************** 
^** ******* *♦**♦*) 
ftKlTcLo»20U0) C IT I f LEL I ) » I =1 » 00 ) 

CUOO FORMATLcri *fdX2h* d0Al»2H *47X1H*) 
mKIT&Lo’ 1999) 

#RITtLt»»20uS)IjT #ItP » IFM #IF0PT » IFTYP1 > IFTYP2» IPLPOSr 
1 DMOoLtiSOPT 1 1ST »KSR »nCIDT »NC1 »NF1 »NF2 » 

c hkIuI *ImK1 »nAVl #CIZPRTCI)»1=1»6) 

*uUS FomMhTC iriU » loX96HXuT IEP IFM IFOPT IFTYP1 IFTYP2 IPLP0S IPM00E ISO 




bAP PKOGRmM LISTING 

a»' 1 iol ..SR UClol nOj NFl NF2 .4KIU T NM NAVE IZPRT/18XU * 2X12 * 3X1 1 ► 
dn Ail* oaa 1 *oX lx* uXH*oa 11*5X11*4X11*3X11*3X1 3*2X13*2X11* 3X11* 2X13* 
3*:X *3* a:aa3*4aI a ♦ 1X1* » IX II » lXll * 1X11 » 1X11 ) 

«KiTi_L o*dGio) Ai*iAX*i>tTA*L)ELT *FSTaRT*FC 
20Xo MJi\t*.AlL jut »OMAi'iAX =»Llb.B»lHr4XfaHULTA =*Llbt8* 1H*4X6HDELT =*E16.8* 
l*h»4Aoiu bl#xi<I = *Elu»8* lli*4X4HFL = *t.lb»6) 

*Fi XFor' i »Lu*u) oO 'l 0 13 
lF'c IF OP I .Lb. 2) 00 To 9 
IKilYPx.wf .3) 00 10 fa 
KtAOlib* lUUb) CF1L I ) * i=l *NF1 ) 
iJub FOKiMaTLoFIu.u) 

i\Laol b* j.uOb) L ToAGxL K) » t=l »NF1 ) 
aLaLil t>» iuob) L TPHlLL) »L-1 hmF 1 ) 

.,MTuLu*2iUl)Ch'il 1) * TMAOlt i) *TPllU *) * 1 = 1 *NF1) 
dUix FORh.hUIH »43A2hFUlXbHTtJ|AGlUX4HTPHl/C35X3E15.8) ) 
uu 1 u 9 
c lOhI iiiUu 

hUAUl o* o.uLib) nOKDl *FCUT 1 »KIP1 
i.Kx 1 clo*2u1u) <0KD1*FCUT1»KIP1 
dOXG J-OKF.aTi. ih »47XahNOKul5X5HFCUTlllX4iiRIPl/50XU»2E16.8) 

9 wOuT 
Xi cOul i.4uu 

xFL IF I YF2.._ T#3) GO TO 12 
XFcIFoPI .t.tt.1) GO TO 13 
kuml-L b» xOUb) L F2L 1 ) 1 1=1 1 NF2) 

Rt*«uLb* XOOb) L li v >AGdLK) »K=1*uF2) 

kEaul b* luiib) L Tph2LL) »L=1»NF2) 

riKxlEL 1 ) t T toAOdt I ) » TPh2C i ) *I=l*hF2) 

2oXd Oi\n*\Tl in 1 4 JAdilPZxl Abh i r<|AG211X4hTpH2/L 35X3E15. 8).) 

GO TO 13 
Id wOnl ImUl. 

At AOt b t i uOb ) U 0 A 02 1 F CEN 2 1 FCUT 2 » RiPd 

10 Ou rOnt'iHlL lb* jFIo » u ) 

WkxTttb*201S) NuRD2»FCEu2»FCUT2*KlP2 
201b FCtthttHin » 35XbhNOi\ Jd * 6XSHF CEn 21 lXbHFCUT 2 12X4HR IP2/36X 1 5 » 3E 16 « 8 ) 

13 CUi* I 1HUL 

lFLlTTPx.Eb.2) oO TO C 14* lb* 14* 16* 14* 16) »N0RU1 

14 cOnI HiUu 

iFLXFT Tt d*i-G»2) GO TO L 17* 16* 17* 16* 17* 16) »N0RD2 
\su TO x7 




I 
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lb COM iwuL 

•VKlTELbr I99d) 
wkilELbrl999) 
rttof 204b) 

204b FOKMATC1H r4bX40HEV£N OROER CHEbYChEV FILTER NOT ALLOWED.) 
itRiTEC ur 1999) 

STOP 

1 7 COM i nut 

iftRITECbr 1999) 
lFLlSOPT.EU.O) GO TO 20 
GO TO L br 10 r 15) r ISOPT 
5 COM iUUt 

KEADl 5r 1010 ) m 
1010 i-OM'iAVL 12) 

hEAuLbriOlb) CALI) »I-lrM) 
kLaOLS* 101b) Lp'L I ) » i=l i k) 
ttEAuL br 101b) L ThETAL 1 ) r 1=1 rM) 
lulS P0 KMm1LoF1u«0) 

>.R a Tc.Co»2U01) LACL) »FCL) » THETAL i.) »l»L=1»M) 

2001 rOKMATC ill r48X4NACh) 12X4HPLM) llX8HlhETAC M) 3XlHM/C42X3E16.8r 13) ) 
*»KaT tC or 1999) 

oO TO 20 

10 COMINUE 

KEhLL 5r 1020 ) TTpCTi r 1=1 » 5J 
KEaLL br 1020) C AL I )j 1=1 r 4) 

1020 FOkMa fL bF10 • 0 ) 

CRITLl b»2002) CTPCD r ACL) »L=lr4) rTpC 5) 

2002 POkMATL lrl rb3X2hm5XlHA/C56X2Elb.0> ) 
wKITfcCor 1999) 

GO To'cO 

lb LOnTlNUt __ _ 

KLA6t5rl025) TPCl) r ACl) 

102b FOKMATL2F1G.O) 

MKITcCor20f 2) TPCl) r AU) 

IPll) = TPLD/2.0 

xRITtCbr 1999) 

20 LOMINUE __ 

GO TO Ci00»200‘r306*400r500r600»700r800»900rl000) rlEP 

xuO CON) 1N0L 

LOcmT = lOo 




I 


SAP PROGRAM LISTING i 

i.itili_Lo*20tU) toCM * luT * IS1 * NRIOT * uC Iul \ 

20c.u t ciM'irt fLinOfHlXooH************************************************* 
i*/HLA2..**4oX2ri**/4 t :x2H**5X30HPROGRMM mAS PASSEO ENTRY POINT* 15* 1H. 
£.DAtn**/‘t2X2li**ox2t>n l)ATA IS CURRENTLY ON UnIT* 12* 1H. 9X2H**/42X2h** 
wl»X 3c-rt jPARE TAPc. IS CURRENTLY ON UNIT * 12* 1H.6X2H**/42X3H** 10H 
4 Tr-iLKt. Arte* itflAttiirtOAS AND * 14 * 1 Xl9riCOLUrti«iS OF DATA. **) 
c a i»l 1*1 ) 

lTi>SwX — M T I Me 

vmLL I h-.c-Kl iTi.-itl * , “iTlrtE#MlN»S£CS) 

«»Ki1 C.L o * 140 j ) .-liN * otCb 1 

l4o3 rvjitr-.A'lL ill * 4lX19ti** EXECUTION TIME =» 13* 1 XcHMInUTES* 2XF6.3* 1X7HSEC 
iui<ub» on **/42X<ln**4oX2H**/42Xb0h********************************** 

t **** ***4-******** ) 

il-Liaupi .Eu.O) GO 10 l5u 

eMuL IgwENL 101 i.^RXu'i *NCIDT *UELT * ISOPT * IFM»a»F*THETA»M»TP) 
iFLi^pKlLl) .t»*.l) Cmi_L I WR1TEL ID! * nCIuT » nRIDT) 
oO Id too 
i bo loM itiUL 

C/.LL ISmRL i J1 * 1ST » hK IDT * NC IDT » OElT * NR1 * NCI » KSR ) 

A = FluMLaSr) 

Ottl = utLl/X 
iTt.r*.p = luf 
1UT - IM 
1ST — i 1 EMP 

iFl i2Pr\ *il 1 f • £<*• 1 ) CAlL TtfRlTEC 1 UT»nCI0T*NRIDT) 

2UI) COM ii'iUu. 

LvLM I" = 20 0 

«KITt.Lo*2U2u> eCCA f » luT * 1ST * WRILT * NCIOT 
ChLL CLOCKL MT ifiE) 

CALL 1 Ai'iERL 1 TiMLX* wTlmE » MIN* SECS) 
iTiWeX - XiTlME 
nRiltlo* 14u3) MiN*^ECS 

c*** uc. f lrmi 4c. MAxiNjh Signal component 

lUM 1-iUT 
iUi»T2 = iST 
NR— Nrtiw I 

uc^NCiur 
bMAX-0 • 0 
uO £. d=l»NC 

call i«iAoTAPC II *i\R* aR* IU.kTI) 
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CmlL l-.Ho I AkL A A » IK » mA » IUn T1 ) 
l-o 1 A— A i uiK 
a-mooL t\><,L A ) ) 

X xFt A «0 I • bPAA) Si»iAX=X 
c lOo I iuUu 

vMuL .'IHL f APC IH MvK t Ml< t IOi\Tl ) 

C*** i"ur\HOi\t‘i uOKnAlIZhT ION 
ai'irtA- A • U/Si'iHX 

UKJ H I wC 

Call xaoImPL ii>i,K»nKf IUinTI) 

CmLL Mr’L 1 a »kK» A A » IU.-Tl ) 

uO b J i«l i i<k 

AM. A )=o. ,«X*At<l. i ) 

Vo cO./Ia.muc- 

l All t«iMo T APL icf uRi HP f XUi iT2 ) 
wAlL MAoTAPL ic»r.K» aI t 1 UnT 2 ) 

A lli\T ihJu 

U»kL I*i«v I aPl l*t»urt t aR» aUuTI ) 

CmlL t-.rtol APL XL»;„K»aR> IUnTl) 

A l Ll*iP= A~< l 
A L> I = A b i 

Abl = ITli'P 

lAll TTRAubL i^KaoT »»<l 1 uT » IUT * IS N 

ITlMP = in 

au) - Abl 
AbT = aTLHP 
Altf-tP = uCAuT 
mCalI = UKal/T 
uttluT = UlmR 
30u CGuTiLUL 

lCCmi = Ouu 

i.KilulurLOLO) lCCa l * IDT 1 1ST * UrtlUT * WCIDT 

CALL CuOCKLi'iTIi>ic.) 

CALL 1 a: ilKL AT XMLX » ,iTlME»MlN»SLCS) 
aIa^lX =■ l*.l Ihl 
t»l<ilLLu*lHu o) i41i4rbt.CS 

aKL1I-oPT.lQ. 1).OR.I:IFOPT.EQ.3>) GO TO 310 
iFl AFi v i*i«£* 1 ) Co TO ouQ 

OAU lOuTaiJUl 

call LH'lLAR»Ai»WRM/IrIuT»IST »NRIOT r hCIOTr II ) 



SAP PROGKA.-I LISTING 


lll>.P - iu'l 
iui - i -j! 
iSi - i i tWP 
iTLfrP - utiuT 
• Xiul = i.KiUT 
..rUtT - iTt.-iP 

itLltPrtiLc) .£o.D gaLl TWR1TEL ICT tUClDT » NR IDT) 

HuU COi.l An Jt. 

uoCAj - GUo 

.>Ki I&.lo*2U2o) uCoAf • iDTi 1ST H .RiL»T» .TCIi/T 
CMt-L. Cl_OCM.Mri.vlt) 

GMtu fi: Lkl IT XMl./ * i'-.TIi'iE » MlN » SECS) 
lliN&A - i*. llrtE 
•<Kl1cL o* l4u.J) MiN»SECS 
C *** * ni I aml. i Zt. i>AlA FOi< FILTER 

uF«tu = l.u/LGEuT *rLuATL(4RluT*NCI0T) ) 
iFLL iFuPT.c.u.i) .OK. LiFOPT. EG. 3) ) Go To HOB 
iFiu - H 
GO To 4lb 
hoc CGi^liNU.- 

iFLlF.vi.tu.i) GO To 412 
ifiL = iFTYPi 
GO TO Alb 
Gi2 COuliiioE 

iFit = iFTYPi + 4 
Hib COi.l A..ot- 

FCwNx - u.U 

A* = *1.0 

CMUU F xtTtKt iOT r ISf ».\C1 i#T»NRIil)T »DFKEG» X9» IFIL»TMAG1»TPH1»F1»NF1 
* 1 1 JOKui » FCEi .1 # FCU r 1 rrtIPl) 
iTcwP - loT 
10T = isT 
1ST = i TEMP 

lFLi..PKTC3) .c.w.1) CAuL TWRITECIOT».iClOT»NRIDT) 
so 0 cOuTIngE 
uOCAT = Sou 

*hiTtLu»2o20) toCAT* 10T» 1ST »NRIDT»NCIDT 
CmLL CuoCKLMTIiml) 

CAgL fiisLKt ITIMEA»MTIME»MIN»SECS) 
iTlNtX " KTitfE 
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»hl'(LLo>14UO) hUN’SLvS 
InOuL = i i 

CALL LFKTtMKr Ai » rfR» wl » IDT» 1ST » NR IDT »NC IDT* MODE) 

iTLhP = IDT 

IDT = 1ST 

ibl = xTEMP 

ITlMP = UCJiuT 

ixCiL/I - NRluT 

NR IDT = ITLMP 

1FL 1ZPK1C4) . £u. 1) CALL T*RITEUUT»NCIDT*NR1DT) 
ouU CONTINUE. 

LOCAT = oOU 

/.KiTL£o*20c.O) coCAT * IDT * 1ST *NkIDT * WCIDT 

Call CluCalmT I.*i£) 

call f ii* thl IT » 1*1 T a ME » Ml SECS) 

1 1 lt*ic.A = i*i I ItoE 

tokllELo*i4ud) MlNrSECS __ 

CAuL KolA IDT » 1ST *NRIDT tivlCIuT »BET a) 

1Tc.I*iP=ilT 

1L1=IS1 

iST-ITti'iP 

70 U COwTlNGt 

lOCAT =700 

hklTLL t>*2UL0) LDCAT#IOt»iSt»NRiDT»NCiDf 
CALL CLOCKCMTImL) 

CALL T X ttEKl I T I ML X * i*.T I ME » M I N * SEC S ) 

1 TIMEX = MTlML 
WklTELb* 1403) MINrSECS 
mOuE =1 _ 

CALL EFFTC Ak»AI »WR» toll* IDT * 1ST #NRlDT»NCli3T » MODE) 
ITlMP = IDT _ 

IDT = 1ST 
i T = 1TEMP 
llfcMP = NCiOT 

.vtCiuf = WRXDT _ _ 

UKlDT = XTtMP 

iFl 1ZPRTC5) •£U*i) CALL TwRITEt IDT »NCIDT»NRIDT) 
dOU COhTINUc. 

lOCAT = dUO __ 

nK1TLLo*LG£u) LOCAT* IDT » 1ST »NRIDT»i)ClDT 
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LhLL Ci_uCNL i*i T l He. ) 

CAuL T ii'itKL i T XMlX t i*l T IrtEiMlN’SECS) 
i I xMla — K f L»ic. 

..RxlELof 14ui) '.1 W » bc.CS 

IK xFUPT *L I ^0 TO 900 

ut-tiuvji = 1 «/l FLGATL i*KI0T*NCXD T ) *DELT ) 

if xL — xFT f P2 

cmuL FxcTEia I^T’IH »uCIiJTHMRIDT» 0FREQ»FC»IFIL»TMAG2»TPH2»F2»NF2 
* » WOUu£ » FCEu2 » F Cb T2 » K 1 P2 ) 
x Tti'tp — iL»r 
io f - 1ST 
xsi = xltMP 

xFl IZPrtlLo) CALL TrtRITEC IDT’NCiDT’NRIDT) 

■i i/ U LOrJ I INuc. 

uOcAf = 90u 

•<KXTELb»2u20) uoC A 1 1 iuT 1 1ST ’NRIDT ’NCIDT 

cALL CLOCKlMf iMfc) 

cacl l IwERt i.riMc.A.»i'iTlME’MlN»SECS) 

IlIMbX = MTlME 

nKxTEL b» i4o3) MiN’SECS 

CAuL TTkAnSCNKXuT mCIDT » IDT » 1ST ) 

HEMP = IDT 

xul - xaT _ 

1ST = HEMP 
xTt.MP = NCI JT 
uClLT - NRIDT 
. nRIDT - I TEMP 
luOU COnTInUc. 

lCCmT - iOUO _ 

»KITlI;o»2020) LGCAT»TdT~» 1ST » NRIDT HjSCIDT 

CALL CLOCKLMriME) / ...... 

C ALL T1 MERL X T I ME X » MT I ME » M I N » SECS/ 

lliMwX = MTIME / 

wRITEth’ 1400) MIN ’SECS / 

L/Fklu — I./luEcT^FLOATLimRIDT/nCIDT) ) 
a = CFS1 R.l ) / uFReQ / 

nstakt = X ♦ 1 / _ 

uLlF = uFRE'gi / 

CALL PlOTL IDT » 1ST »r«RIL f »NCI0 T’ IPmOJE’ NAVE»FC»DELF»FSTRT 

1» AmmA» XPLPOS) 
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oRiTELo»1999) 
wRITliL of<iOcb) h!>iAX 
tLvtb POkKim Ti. xh » orlAi-ifiX = ft,15.8) 
r.KlTdLc»1999) 

iK IPMooE.LG. 2) GO TO 1050 
*«RlTcL o» 20.}$) 

iiuib^POKNATLlh ►HOAbuNRON TERMINATED. PLOT DATA TAPE GENERATED ON UNIT 

oo ro iioo 

lObU COhHhuc 

r.KlTdL or 2040 ) 

2GH0 POkMaTLIH »43X44HRUH TERMINATED. NO PLOT DATA TAPE GENERATED.) 
lioO cuuTIijoc. 

wkildCbrl999) 

CALL CuOCKL ITIimu2) 

call TlistRL lTiMEir iTIME 2»M1N»SECS) 

aMTcC or 14ub) Mii.rSe.CS 

SlOP 

Lf*U 


SAP program listing 


Cd* 

CD 

CL 

CL 

CL 

CL 

CL 

CL 

CL 

CL 

CL 

Cd 

CD 

CL 

CL 

CL 

CL 

CL 

CL 

CL 

CD 

CD 

CL 

CL 

CL 

CL 

CL 

CL 

CL 

CD 

CL 

CL 

CL 

CD 

CL 

CD 

CL 

CL 


gUBRLU ( iNL PLO U IU., ( l # IUNT* » NR » NC » i PMODE t NAVE t FC t DELF » FSTRT 

1* Ai'iAa* IPLPOS) 

********************************** 

PROGRAMMER ANL LATE 

KiCHAKu C • THOi'iAS 
TrtW SYSTEMS 
JUNE 1*72 

rLrtPoSE 

GENERA fES THE LATA TAPE REQUIREO TO DRIVE THE EAI PLOTTER 
ALSO PRINTS The. DATA TO BE PLOTTED. 

USAGE 

CAuL PluTl I UNT 1 » 1UNT2 1 NR » NC » I PMODE » NAVE » FC » DELF # FSTRT » 
*MAX* IPLPOS) 

uESCKIPTIO.m OF PARhMETERS 
luPUl 

CALLIng SEgUEnCE 

iuuti - input tape unit 

IOnT2 - SPARE TAPE UNIT 

NR - NUMBER OF ROWS IN INPUT ARRAYS 

NC - NUMBER OF COLUMNS IN INPUT ARRAYS 

NAVE - NUMBER OF INPUT POINTS TO BE AVERAGED IN OUTPUT 
FC - CENTER FREQUENCY OF PLOT 
“ DElF - INCREMENT BETWEEN INPUT POINTS! FREQUENCY) 

FSTRT_ - STARTING FREQUENCY OF PLOT_ ___ 

IPMOLE .EQ.O PRINT POWER SPECTRUM AND GENERATE DATA 

.EQ.2 ONLY PRINT POwER SPECTRUM 

•EQ.l ONLY GENERATE DATA TAPE 
AmaX - MAXIMUM POWER COMPONENT OF INPUT SIGNAL 
COi’iMON 

I TITLE - UP TO 30 CHARACTERS FOR RUN TITLE AND/OR DATE 

CARL 

Tape 

AR - REAL PART OF SIGNAL 

AI - IMAGINARY PART OF SIGNAL . _ . 


2-61 



SAP PROGRAM LISTING 


CL OUTPUT 

Cu CAUL bc.tfUc.iXE 

Lb COMMON 

CU PRINT 

Cb P POwEK SPECTRUM ARRAY » DEPENDENT PLOT VARIABLE 

Cu F FKEUUtNCY ARRAY » INDEPENDENT PLOT VARIABLE 

CD TAPE 

Cu 

Cb REmAkkS ANl» restrictions 

Cb SUbKuUTiNtb REOUIRlD 

CU MA6T*P 

Cb 

Cb* 1 ********************************* 

UlMLNSiON MRCiSc) »Alt2bte) »Pt2U46> »FC2048) 
uIMtNSiON 1 TiTutC 3u ) 

C0MNON/STun6£/aK» A i »P»F 
COMMON/ 1 ITl/ITITLE 
11 = 1 
12 = 2 
a4 = 4 

ib = 5 _ 

IS12 s 4«0 

CALL MAoTAPL 14» I4r AKf IU nTI) _ 

CAlL KaGTAPL 14» 14* AR~» IUnT2) ~ 

IFtAMAX.GTiO.U) 60 TO b 

UO 4 J~i * |«C 

CALL MAuTAPI I1*wR*AR* IUNT1) 

CALL MmuTAPL II hmK* Al * IUuTl ) 

uO O 1- J. * Nr\ _ 

X = ARCi)**2 + AlVl)**2 

1FLX*6V . AMAX) A,*tAX = X „ 

i CONTINUE 
4 CONTINUE 

CALL MAUTAPI. 14* 14* AK* lUiiTl) 

b CONTInUl 

UPS = 4ttli 

IFL IPmuuE»£u»2) nPS = NR*NC/2 _ 

iFLFSTKT) o»7»7 
b CONTInUl 
1FC = 0 
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uu To 9 
7 COt.T 1NUL 

IftFSTKi.Gt.FC) GO Tv 8 
1FC = -a 
oO TO 9 
b CONTINUE 
Aft = 1 
9 lOnTAiwl. 

IFL IFC.i.t.O) GO TO 20 
*PaSS = 1 
rtPi = NPS/2 

FkaX = v£lF*FL0ATL nR*NC/2-1 ) 

FLAST = UElF*FlOATCnAVE*CNPS- 1) ) 
iFLftASl.Lt.FMAX) GO TO 30 
NAVE = LLNK*NC/2-l)/Ci'lPS-l) ) 
ttHITtt or 2U0U ) nAVE 

2000 FOKMAU1H r 25HuAVt HAS BEEN ADJUSTED TO* 14) 

GO TO 3u 
20 CONTINUE 

FLAST = FSTKT ♦ DEtF *FLOATC NA VE*C NPS-1 ) ) 
lFUFC.EG.l) GO TO 22 
FMAX = FC 

*0 TO 24 

22 COnT lNUt 

FMAX = FC ♦ 0ELF*FL0ATCNR*NC/2-l) _ 

24 CONTINUE 

if L FLAST. LE.FM mX) <jO TO 30 

UPS =t INTL L FmAX-FSTKT ) / wtLF ) 41 ) /NAVE 

*R;lEL6r20u5) nPS 

2j 05 F0RMATC1H *45HNUMd£R OF POINTS PLOTTED HAS BEEN AOJUSTEO TO* 14) 
iFLNFS.GE.4dii) GO TO 30 _ . 

Iti Abt = NPS + 1 
UO 2b I=IBASt»480 
PCI) = -70.0 

2b CONTINUE _ . _ 

30 CONTINUE 

EPS ** 1 . E-7 _ . . . . . . .. _ _ 

SOM - 0. 

nSuM =0 

A = i 


| REPRODUCIBILITY OF THE ORIGINAL PAGE. IS POOR 
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PMIN - ”70 • 

AAVL = uAVt 
it-LlKC) 5502*04 
Sc CONl INUt 
u = u 
nT=1 

l-S l>\T - FC 
vu To oo 

54 COM iiMOL. 

uSfMKT = 1 NTLlFSTRT”FC+E.PS) /DELF ) ♦ 1 

u = NSlART/iv»R 

isl = NoTAKT - h*NR 

wO I sJ Jc 
JO AiVJu 

t UUN - KC - DtuF*FL0*TCNR*hC/2) 
i.SImkT = lNTCLFbTKI-FGUM+EPS)/DELF) + 1 
14 = imSi'AKT/NR 
nT = USlAKl - N*NK 
N = u + i>*K /2 
id couiXhUt. 

40 CONTINUE . . 1 

IFLN»LU*U) GO To 5U 

uO 45 l—l>u _ 

CALL MAGTAPCU»NR*MR*lUr4Tl) 

CALL MhGTAPLI1»nK*AI»IUNT1' 

45 COnTINOL 
5u COuT iNUt 

CALL MAGI ARC ll»NK» AK» IU nTI) 

CALL MhGTaPL II »NR*AI* IUNTI) . _ 

uO b5 1-1 »NK 

ARLI) = ARE. 1 ) **2 ♦ Alt I )**2 

55 CunTiUOl 
f-U COnT InOl 

1FL N51M*EGi*0) SUM - 0* 

SUM = 5UM + AKCkT) 

AT = K1 ♦ 1 
i«SUK = hSUm ♦ 1 
IFLUdUM*LT*NAV£) GO TO 70 
KH - K 

IFUIFC.EQ.OI.ANU.CIPASS.EQ.I)) kp = KP ♦ ‘nps 
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«R»/ — SuM/L Ai‘1«A»XAVL) 

aFLAKG.ll.U.O) AKG = l.Ot-08 

FlKP) = 10.0 *AllG1uLAKG) 

iFtPLAP) .LT.PMln) PLkP> = PM IN 

FlKP) = FSTRT ♦ DEwF*FLOATLCK-1>*NaVE) 

N = K + 1 

,«SUf*. s o 

Ii-lh.GT.NPS) GO TO oU 

7U COmTINuL __ 

1FLKT.LL.NK) GO TO ou 
kT = 1 
bO To DU 

oU COimI ahUl . _ 

IFLlFC.Nt.U) GO TO 90 
iPMbi = IPASS + 1 
AFLlPASS.GT.a) GQ TO 90 

SOn — U.O 

NSUM : 0 
K - 1 

NSTART = NK«riC + 1 - nPS*NAVE 

n - NSTAKT/NR . — 

KT = NSTART - N*NR 

FSTRT = FC - OELF*FLOATLNPS*NAVE) . _ 

CALL KA&TAPC I*t»I4»AR» iUwTl) 

GO T* 40 .... 

9C COUT AhUl 

IFL IPnODE.lo. 1 ) GO TO 100 . 

irtIFC.fcQ.0) NPS5H8J 

nPSI = (mPS/4 

hPS2 = 2*NPS1 

NPS3 = 3*NPS1 . - 

ftRXlcCbt 1999) 

1999 FOKnATC 1H » I32rt»*********«t**************************************** 

X#*****»*********»************************************************* 

£***************) _ _ 

rR11£Co* 2010) C Ft K ) »PCK) »FLK+NPS1) . PtK+NPSl) #FTK+NPS2) »PCK+NPS2> # 

1 FLk+NPS3)»PLM-NPS3)»K=1»NPS1> .... . 

2010 FOKMATL lril » 5X4r1Ft K ) 12X4HPC K) 13X4HFC K ) 12X4HPC K ) 13X4HFC K ) 12X4HPC K ) 13 
1X4HFL K ) 12X4MPC K ) /C 2E16 . 6 » 1X2E16 • 8 » 1X2E16 • 8 # 1X2E16 . 6 * ) _. 

»RiTEf6»1999) 
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iFLlPMOl/t.cU.*) GO To 110 
100 CCNTlNut 

HI) = 1PLPJS 
uO 105 l=lrJ0 
FLi+1) = ITITLLl 1 ) 

105 cGNTlwut 

FL 02) = DEuF 
I-LjO) = FC 
Ft 04) = FSTRT 
FLOS) = NPS 
FLO b) = IFC 
lOo - Oo 

CALL Mm*TAPCI4»14»P#I2) 

CALL mAGTAPL Ic.# I0b»F # 12) 

CALL MGTAPC 12# l512»P# 12) 

CALL MAGTAPC I5#14»P» 12) 

Call hAGTAPC 14 # i4#P» 12) 

110 COiNlXNOL 
Rtl URi\ 

ENw 
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C 




><> 


>/’* I 


* 

*A . • 

>V 
t-’ ••' - 

• >r! . 


CD* 

Cu 

CD 

CD 

CD 

CL 

CD 

CD 

CL 

CL 

CD 

CD 

CL 

CL 

CL 

CD 

CD 

CD 

CD 

CD 

CO 

CL 

CD 

CD 

CD 

CL 

CD 

CD 

CD 

CO 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 


SAP PROGRAM LISTING 

SC’OKDDl iiMt iii IReVLr.lXX.KP) 

L - i_XX 

********************************** 

PROGRAMMER ANu DATL 

OAt/ lo. m. UcTCHMENUY 
Tk« SYSYti’O 
MaY 1*70 


POrtPobt 


Tills subroutine solves the bit reversal problem » I.E. 

GiVe* R AND L aHEKE 

K = iu + I1*C2**1> + ... ♦ lL*t 2**L) 

Find nP ahere 

rP = I0*C2**L) ♦ Il*C2**CL~l) # + ... ♦ IL 
IN THE ABOVE I0» U» ... » IL ARE ZERO OR ONE 


cSaGe 


call oITREVLml.kp) 

DESCRIPTION Or PARAMETERS 
XNPol 

calling sequence 

k - integer 

l - defines number of binary digits In k 
Common _ . 

NuNE 

Card _ ..... 

none 

Tape 

none 

OUTPUT " .. 

calling sequence 

KP - INTEGER RESULTING FROM THE BIT REVERSAL 
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Co 

COPiMOw 


CO 

NONE 


cu 

ChKO 


cu 

NO.«E 


Cu 

PKiKl 


cu 

NOhl 


cu 

TmFE 


cu 

cu 

NONE 


cu 

RLMAKKS AWU RESTRICT IONS 


cu 

cu 

NONt 


CO 

SUoROOriNtS KLuuIKlJ 


CO 

CO 

NOnE 


CO 

totTHOO 


CO 

TESTS EACH BIT IN K BY 

SUBTRACTION 

CO 

OF TH£ K£SUi_T TO GENERATE A 1 OR A 

CO 

CO 

BINARY OIGIT OF KP 


CO* 

***************** 
lO = 2 

****** 


RP = 0 
R0 = R 

. . . 


01 = 1 
bi> - 2**L 
UP = L+l 
O0 10 i— 1 »LP 
RN = KO - 02 



1FL RN .LT. 0 ) GO TO 5 



RP = RP ♦ 01 
RO = Kn 



b J1 = 2*J1 
10 J2 = J2/2 
rttTUKN 
LNb 
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CD* * 

CD 

Cu 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

Cu 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD ’ 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD* * 


SUBROUTINE ClOCkCJ) 

********************************* 

PROGRAMMER AND DATt 

RICHARD C. THOMAS 
Truk SYSTEMS 
JUNE IS72 

PUKPOSl 

ClOCK IS A DUMMY ROUTINE THAT ALLOCATES CORE STORAGE FOR 
A CS-i ASSEMBLY LANGUAGE ROUTINE NAMED CLOCK. THE CS-1 

program references the computer clock, the computer clock 

GIVES AN ABSOLUTE TIME WHICH HAS NO MEANINGFUL REFERENCE 
TIME dASt. TwO SEPERATE CALLS TO THIS ROUTINE KAY BE USED 
WITH SUBROUTINE TIMER TO COMPUTE THE RUN TIME OF A SEG- 
MENT OF THE PROGRAM. 

OSAGE 

Call cloCrLi) 

DESCRIPTION OF PARAMETERS 


iNPuT 

NONE 

OoTPuT 

Call SEQUENCE 

I - TIME FROM COMPUTER CLOCK. 

REMARKS AND RESTRICTIONS ” " 

THIS ROUTINE ALLOCATES THE NECESSARY CORE REQUIRED FOR 
THE CS-1 PROGRAM WHICH MUST BE LOADED BY PAPER TAPE. SEE 
PROGRAM OPERATING INSTRUCTIONS FOR CORRECT PROCEDURES. 
SUBROUTINES REQUIRED 
NONE ' 

************** * * ***************** 

DIMENSION AC 300) 

aC1)s1.u 

RETURN 
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SUoKoUT iNE EFF lL AR » Al » WK > wl f 1 1 » 10 » NRO*» NCOL » MODE) 

CD* ********************************** 

uO PKOtoRAftMfcR AND BATE 

CO KlCilARo C • THOMAS 

CO Tk.> SYSTEMS 

CO JtiCLrtBtK 1971 

CO 

CO PURPobfc. 

co computes the discrete Fourier transform using the cooley- 

co TUKtY algor iTHM 

co 

CO oSmGE 

CO CALL EFFTL Ak»A1#WR»1«I»II»IO»NRO«» NROv» t NCOL # MODE # NBUF ) 

CO 

CO uLgCkIPTIOn OF PARAMETERS 

CO 

co Input 

CO AH - Rc.Au PORI ION OF FUNCTION TO BE TRANSFORMED. 

CO Ml - IiMAGINmRY PORTION OF FUNCTION TO BE TRANSFORMED. 

CO - Rt.Au PORTION OF EXPONENTIAL TERMS. 

CO wi - imaginary portion of EXPONENTIAL terms. 

co li - input tapl unit. 

CO 10 - OUTPUT tape unit. 

CO NROfc - NUMBER OF ROlfkS IN INPUT ARRAY. 

CO NCOL - NUMBER OF COLUMNS IN INPUT ARRAY. 

CO MOut - .EW.l COMPUTE DIRECT FFT . 

Co .NE.i COMPUTE INVERSE FFT. 

CO 

CD 0U1PUI 

CD AR - RLAL PORTION OF TRANSFORM. 

Co Al - IMAGINARY PORTION OF TRANSFORM. 

Co wR - REAL PORTION OF EXPONENTIAL TERMS. 

CL *1 - IMAGINARY PORTION OF EXPONENTIAL TERMS. 

Co 

CD dUdKOUT *NtS REQUIRED 
CU aiTREV 

CD FFT 

CD* ****************** w *************** 

COMKON/STORGE/ AR 'Ai’WK’Wl 

DIMENSION ARC 266) » AIC256) »rtRC204fi) »ttIC204fl) 
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! 

, id = 2 

14 = 4 
imK = NKu* 
hC = NCOL 

, CC'W = 6.i'ttiid54/CNR*NC) 

' CALL rtACTAPL 14* I4 *mR» II ) 

CAuL tiAbTAPt 14* l4*«K* 10) 

. 0 - i4K 

N = 1 

10 COul IimuL 

lFLO»Lu*4) bO To 20 

. I» " U + 1 

i O = 0/2 

GO 1 U 10 
20 COwTlNUL 

r\ if-LKOLL.EG.i) CON = - CON 

W UO 40 0=1 »NC 

CALL NAGT APt 1 1 * wK * hk * 1 1 ) 

CALL MAGT APt 1 1 * NR * A i * 1 1 ) 

CALL FFTLN*M0LL*0) 
uO 50 l=l»Nrt 

A = CON*FLOATtU-l)*CJ-l)l 

Twk = COSt A) 

Twl = SlNtA) 

ARP = aRCI) 

A1P = AILI) 

ARC!) = ARP*TwR-AlP*TwI 

AILI) = aRP*Tw1+AIP»TwR 

30 CONTINUE. 

CALL MAGTAPtI2*NR»AR*I0) . . 

CALL MAGTAPt I2»NR»Al» 10) 

40 COuTlNUc. 

CALL M AbT APt 14*14* AR» II) 

CALL HA6TAPtI4»l4*AR»I0) __ 

CalL TTKANSt NR » NC * 10 r 1 1 ) 

L = NC 

U - 1 

SO COuTINUl 

IFLL.fcu.4) GO TO 60 


c 


2-72 


SAP PROGRAM LISTING 


l<i = IM ♦ I 
U = U/2 

wO lu SO 
60 CONTINUE. 

LG 7U J-I»wR 
CAuL MAGTAPC H»dC» AR» II) 
CALL MAGT API 1 1 r UC * m I » 1 1 ) 
CALL FFTCN»MUUt» J) 

Caul MAGTAPt I2 »nC» ak» 10) 
CALL MAufAPL Ik'NCrAlrlO) 
70 CONl I Nut 

CAUL MAGTAPl 1A» I4 »mK» II ) 
CAuL riAOTAPC I4» 1H»AR» 10) 
kLTOkN 

Lhu 
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bUoKLUl iNt FFTL e.P » MOLE » NCMP ) 

CL* ********************************** 


C u 
CL 
CL 
CL 
CL 
CL 
CL 
CL 
Cl 
CL 
CL 
CL 
CL 
C U 
CL 
CL 
CL 
CL 
CL 
CL 
CL 
CL 
CL 
CL 
CD 
CL 
CL 
CO 
CL 
CL 
CO 
CL 
CL 
CL 
CO 
CL 
CL 
CU 


PROGRAMMER AND LATC 

DAVID i».» uETCHmENLY 

TK« systems 

MAT 1970 


PLKPOSl 

CuiMPUTES Tnt DISCRETE FOURIER TRANSFORM USING THE COOLEY- 
TDHEY mLGuKITHM 

uSAGl 

CALL FR7CAK»AI»«'R»WI»M»M0UE»NCW) 
utSLKXPTiON OF PARAMETERS 
InPuI 

CalL I NO SEQUENCE 

AK - REAL PORTION OF FUNCTION TO BE TRANSFORMED 
AI - IMAGINARY PORTION OF FUNCTION TO BE TRANSFoRMEO 
»K - REAL PORTION OF EXPONENTIAL TERMS 
Ml - IMAGINARY PORTION OF EXPONENTIAL TERMS 
K - DEFINES DIMENSION OF AR» AI» MR’ AND Wl ALL ARE OF 
DIMENSION 2**CM*1) 

M Out - .EO.l COMPUTE DIRECT TRANSFORM 
.NE.l COMPUTE INVERSE TRANSFORM 
NC* - .£0.1 COMPUTE mR AnD mI ON ENTRY TO FFT 

• NE.l WK AND Ml ARE INPUT * USUALLY COMPUTED ON A 
PREVIOUS ENTRY 
COMMON _ 

NOME 

CARD 

NOME 

TAPE 

NOME 

OUTPUT 

CALLING SEQUENCE 

AH - REAL PORTION OF TRANSFORM 
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CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

Cu 

CD 

CD 

CD 

CD 

CD 

Cu 

CD 

CD 

CD 

CD 

Cu 

CD 

CD* * 


C 

C 

C 

C 


Ml - IMAGINARY portion of transform 
WR - REAL PORTION OF EXPOuENTIAL TERMS 
W 1 - IMAGINARY PORTION OF EXPONENTIAL TERMS 
CUi-iMOi* 

NONE 

CAhD 

NONE - 
Hr* NT 

NONE 

TmHE 

NONE 

REMARKS AND htSiKICTiONS 

1. ODTHUT AR AND AI ARE OVER WRITTEN ON INPUT AR AND II 

2. THE OUTPUT wR AND W1 ARE STORED IN BIT REVERSED ORDER 

3. AR AND AI ARE IN NATURAL ORDER 

4. AR#hI#aR» AND Wl ARE EXTERNALLY DIMENSIONED 

SUuRwUl luES REQUIRED 
bITRtV 


METHOD 

REFERENCE - D.M.DETCHMENDY# ENGINEERING DESCRIPTION OF A 
FAST FOURIER TRANSFORM ALGORITHM# TRW IOC 
5522.7-70-113# JULY 1970. 

********************************* 

uImEnSION ARC 25b) »AIC256) #WKC2046) #WlC2U4fi) 

COmMOn/STORGE/ AR # Ai » nR# a I 
M = MP 
NCW = NCwP 
N = 2**LM+1) 

PI = 3.1415927 
ITwOri = 2**M 
CON = Pl/ITinOM 
MK1 = M-l 


generate cosines and sines in bit reversed order# if 

NCW = 1 
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1FL uLm • Nil* 1 ) Gu TO i.2 
uG io i=l*ITrtOrt 
INI = 1-1 


CALL DiTKEVt Itoi *MMJ »KP) 
LP = KP 

X = CON*FLOAUlP) 

«RLi) = COSlX) 

10 MIL I ) = SlNtX) 


12 uGnTINUc. 


SCALE AND FORM CONJUGATE* IF MODE = 1 


MODE .NE. 1 ) GO TO 22 

UFn = i.o/Floatln) 

uO 20 i=l*u 
AM. I) - DFi\l*Aktl> 

20 Altl) r. -UFN*AIU) 

22 CONTINUE 


BEGIN FFT ALGORITHM 

THE INDICES KEEP TRACK OF 

11 - THE A BEING CALCULATED IS A SUB II 

12 - THE GROUP OF A SUB II BEING CALCULA.ED 

IS 12 

13 - THE INDEX ON THE A PAIR BEING CALCULATED 


mPl s rt+1 

DO 100 11 = 1»MP1 

IIP = 2**CI1-1, 

12P = 2**LMPl-il) 
1LXTP * 2* I2K 
DO BO 12 = 1*4 IP 
1EXT = CI2-1)*JU:XTP 
*3R = wrtC 12) 

*31 = »• XC 12) 

DO BO 13*1 * 12P 
1N0 s 13 ♦ I EXT 
INI = 1N0 ♦ 12P 
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AOuK = AKCiuO) 

AOul = AILIuU) 

AOIK = £ kL I hl ) 

AOII = AiClUl) 

CFR = *JK*AUlK - *31*A01I 
CFI = wJI*A01R + W3R*A01I 
ARtlNO) = HOOrt ♦ CFR 

AlClNO) = AOOl + CFI 

ARC lNl) = AGOR - CFR 

bO aICIimI) = AOOl - CFI 

60 COimT INuc. 

100 CONTINUE. 

INTERCHANGE The A-S USING BIT REVERSAL 


mM = M 

UO 120 A = l*w 
AMI = N-l 

CAeL b 1 7 RE vX AMI * MM * KP ) 

IF CAP - KM1) 120*120*114 
114 AOok = ARCA) 

AOOl = AILK) 

AOlR = ARCAP-H) 

AOl 1 = AXCaP+1) 
mKLK) = AOlR 
AICK) = AOII 
aKIAP+ 1) S AOOK 
AlCKP+1) = AOOl 
120 CONTINUE 


FORM CONJUGATE* IF MODE =1 

IFC MOut #NE. 1 ) GO TO 124 
DO 122 1=1 *N 
122 AIL 1 ) = - AICI) __ 

124 CONTINUE 

RETURN 

END 


SAP PROGRAM LISTING 

SUdRoUIINE FILTcRC iTrtPl t ITAP2»NC»NR »DFREQ*FC * IFIL»TMAGrTPH»F#NF t 
1NORDER* FCLh'FCUThUP) 


CD* 


**************** r: ********* 

CO 



CU 

PROGRAMMER and 

OATE 

CO 

KACHARD C. THOMAS 

CO 

TRW SYSTEMS 

CO 

JUNE 

1972 

CO 



CO 

PUKPOSt. 


CO 

performs the filtering operation by MULTIPLICATION of the 

CO 

TRANSFORM OF THE signal BY The FILTER TRANSFER FUNCTION. 

CO 

THE Fll.TER TRANSFER FUNCTIONS CAN BE INPUT BY THE USER OR 

CO 

IT CA 

:« BC SPECIFIED AS ONE OF THE BUILT-IN BUTTERWORTH OR 

CO 

CHE ■’ 

CHEV FILTERS. THIS ROUTINE IS ALSO USED TO 9NT57R1T5 

CO 

THE ^ 

IGNAL» BY THE USE OF THE TRANSFER FUNCTION 1/CJW)» 

CO 

IN THE CASE OF FM SIGNALS. 

CO 



CD 

uSAGE 


CO 

CALL 

FILTER! ITAP1 » ITAP2 » NC » NR » DFREQ » FC * IFIL * TMAG »TPH»F» 

CO 


NF r NORDER » FCE.x # FCUT » R IP ) 

CO 



CD 

DESCRIPTION OF 

PARAMETERS 

CO 



CO 

INPUT 


CO 

CALLING SEQUENCE 

CO 

1TAP1 

- INPUT TAPE UNIT 

CO 

ITAP2 

- OUTPUT TAPE UNIT 

CO 

NC 

- NUMBER OF COLUMNS IN ARRAYS 

CD 

NR 

- NUMBER OP ROMS IN ARRAYS 

CD 

ofreq 

- FREQUENCY INCREMENT 

CD 

FC 

- CENTER FREQUENCY 

CD 

ifil 

- .EQ.l BUTTERWORTH FILTER 

CD 


- *EQ*2 CHEBYCHEV 

CO 


- .EQ.3 INPUT DATA FILTER 

CD 


- .EQ.4 INTEGRATION FILTER 

CD 


- .EQ.5 PRODUCT OF 1 ANO 4 

CO 


- .EQ.6 PROOUCT OF 2 AND 4 

CO 


- .EQ.7 PRODUCT OF 3 ANO 4 

CD 

TMAG 

- MAGNITUDE OF INPUT FILTER TRANSFER FUNCTION 



s . a ' 


* . -V , 


O 


y-* 1 

ji c 








V • :.;v. -v. ' 

r^n.*..-S 


;■< 

r~ 


_jS 


CD 

TPH - i 

CD 

F - 1 

CD 

NF - 1 

CD 

NORD 

. CD 

FCEN - < 

CD 

FCUT - i 

CD 

RIP - l 

CD 

common 

CD 


CD 

AX 

CD 

w:. - i 

CD 

wl - 1 

CD 

CARD 

CD 

TAPE 

CD 


CD 

OUTPUT 

CD 

CALL SEQUENCE 

CD 

COMMON 

CD 

CARD 

CD 

PRINT 

CD 

CD 

tape 

CD 

SUBROUTINES REQUIRi 

CD 

MAGTAP 

CD 

thfn 

CD 

CD . 

NONE 


SAP PROGRAM LISTING 
CORRESPONDING PHaSEEIN RADIANS) 

FREQUENCY! In HERTZ) AT WHICH MAGNITUDE AND PHASE 
DIMENSION OF TMAG»TPH AND F 
ORDER OF THE bUTTERWORTH OR CHEBYCHEV FILTER 
CENTER FREQUENCY FOR BUTTERWORTH OR CHEBYCHEV 
CUTOFF FREQUENCY FOR BUTTERWORTH OR CHEBYCHEV 
RIPPLL^FOR BUTTERWORTH FILTER 

REAL PART OF FUNCTION BEING FILTERED 
IMAGINARY PART OF FUNCTION BEING FILTERED 
DUMMY COMMON VARIABLE 
DUMMY COMMON VARIABLE 


C*«** 


****** 

DIMENSION FC 1 ) f TMAGC 1 ) t TPHE I ) _ 

DIMENSION ARC25G) * AlE25b) » WRE2048) #WIE2048) 

COMMON/STORGE/AR»AI»WR»Wl _ 

11=1 

12=2 _ __ 

14=4 

IFC=0 - 

IFCFC.GT.Q.OI) 1FC=1 
CALL MAGTAPE 14 »NR»AR»ITAP1) 

CALL MAGTAPE X4»nR» AR» ITAP2) 

SET XFdAG AND IOPT . . __ 

IFEIFIL.GT.3) GO TO 1 


****************** 
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iFuAo=x 

iOPT=IFI_ 

GO TO 2 

1 1FLAG=3 
IOPT=IFlL-4 

IFC IFIL.EG.4) 1FLAG—2 

2 COnTAMJE, 

1FU1FIL.EQ.3).0R.CIFIL.EQ.7>) GO TO 4 
uO 1U 1G 

4 DElF=FLO AT t C WC*hR/2 ) -1 > *DFREQ 
lFLFC.uT.0.01) GO TO 6 
FMIN=FC-UELF 
FMAX=FC+uEi-F 
1FC=1 
GO TO 7 

6 FMIN=Q.U 
FMAA=UCLF 
IFC=0 

7 IFLtFMlN.LT.FCl)). OR. CFMAX.GT.FC nF)>> GO To 997 
10 u=l 

3 CONTINUE. 

CALL mAGT APC 1 1 » NR » AR f 1 TAPI ) 

CALL MAGTAPC II »NRr Al' ITAP1) 

NS=1 " ' ' *' - " 

IFC 1FLAG.EQ.2) GO TO 40U 
C*** FILTER aR AND AI FOR FREQUENCIES ABOVE FC 
IUP=NR/2 
00 1U1 1=1»IUP 

FREQ=FLuATC C J-l ) +C 1-1 ) »NC ) ♦DFREO + FC 

CAlL TRFNt FREQ * IOPTVnORDER ~t FCEN » FCUT t R IP » TMAG * F • NF * NS »TR # T 1 1 TPH) 
ARP=ARC 1 ) _ 

AIP=AIC1) 

ARLI)=ARP*TR - A1P*TI 
AICI)=aRP*TI ♦ AIP*TR 
101 CONTINUE 

NS=1 ' '■ ‘ 

IFC IFC.nE.u) GO TO 300 __ 

C*** FILTER AH and AI FOR NEGAtlVE FREQUENCIES C MODULATING SIGNAL) 

IS= NR/ 2 41 
00 201 I=ISrNK 
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iP= NR ♦ UNR/2)+l)-I 

FRfcU=FLOATL NK*.,iC-C o-l ) -l IP-1 ) *NC ) *DFREQ 

CALL 1 r\F Nt FKEw # iOPT » NQRUER » FCEN* FCUT * RIP * TMAG » F » NF * NS * TR » TI * TPH ) 
AkP=AKLIP) 

AiP=MlCiP) 

TI--TI 

ARUP)=mkP*TR -aIP*T1 
A1LIP)=aRP*T1 +aIP*TR 
201 CONTINUE. 

1FIIFLAG.EG.5> GO TO 400 
faO TO 500 

C**» FIuTtR mR AND AI FOR FREQUENCIES BELOW FC C MODULATED SIGNAL) 

500 COimI InOE 
IS= hK/ 2 + 1 
00 5U1 I=IS»NK 
lP=NK + NR/ 2 ♦ 1 - i 

FREQ = F C-FlOATL NR*iNC-C J-l ) -C IP-1 ) *NC ) *DFREQ 

CALL TRFnL FREQ » 1 OP f » NORDER » FCEN » FCUT * R I P * TMAG » F » NF * NS » TR * T I » TPH > 
ARP=aRliP> 

A1P=mIL1P) 

hRL1P)=ARP*TR - A1P*TI 

501 AUIP)=ARP*TI ♦ AIP*TR 

60 TO 500 _ 

C*** INTEGKmIIOn FILTER - MODULATING SIGNAL ONLY 
400 CONTlNUL 
IUP=NR/2 
DO 415 I=1*IUP 

FREQ = FLOATCt J-l)+C 1-1 ) *NC) *DFREQ 
IFCFREQ.LT.1.0E-50) 60 TO 405 
TI=-1./FKEQ 
GO TO 410 
405 CONTINUE 

II = 0 . u 

410 CONTINUE 
ARP=ARC I ) 

AIP-Alt i) 

ARU )=-AIP*Tl 
Alt 1) = ARP«T1 

415 CONTINUE _ „• . 

1S= NR/2 ♦ 1 
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U> 430 I=1S*NR 

FRfcG=FuOATC IMC*NR-C 0-1 > -C 1-1 ) *NC ) *OFREQ 
lFLFREw.LT. 1.0£-30) GO TO 420 
TI= l./FREw 
GO TO 425 
420 CONI INoE 
Tl = O.u 
425 CONTINUE 
ARP=AKC 1 ) 

AlP=AICi) 

AKtIJ=-AIP*TI 
Aid) = AKP*TI 
430 CONTINUE 
C*“** WRITE TAPE 

1FUIFC.NL.QJ .OR.CJ.NE.l) ) GO TO 501 
ARID =0.0 
All!) = 0.0 
501 CONTINUE 
50 0 CONTINUE 

CALL MAGT APL 12 » NR » AR » IT AP2 ) 

CALL MA6TAPC I2»NR » AI # ITAP2) 

J=J*1 

IFCJ.L£.NC)_ GO TO 3 _ 

CALL MAGTAPC 14»NR» ARr ITAP1) 

CALL MAGTAPCI4»hR»AR»IIAP2) . . 

RETURN 

997 *RITECo»998) „ . 

998 FORMAT! 2X » 49H INPUT FILTER FREQUENCY SPREAD INVALID FOR SIGNAL) 

STOP _ 

END 
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CD* * 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD* * 


SUBROUTINE MODC IUNT 1 » IUNT2 » NR » NC # BETA ) 

********************************* 

PROGRAMMER AND DATE 

RICHARD C. THOMAS 
TRw SYSTEMS 

DONE 1972 


PURPOSE 


PERFORMS THE EXPONENTIATION OPERATION USED TO GENERATE 

the modulated signal representajion from the modulating 

SIGNAl 


usage 


- . CALL M0DCIUNT1»IUNT2»NR»NC»BETA) 
DESCRIPTION QF PARAMETERS 


input 

CALLING SEQUENCE 


IUNT1 - INPUT DATA TAPE UNIT _ . 

IUNT2 - OUTPUT DATA TAPE UNIT 

NR - NUMBER OF WORDS/RECORD 

NC ▼ NUMBER OF REAL RECORDS 

BETA_ - MODULATION INDEX# MULTIPLIES ENTIRE TEST SIGNAL 
COMMON 

AR - REAL PART OF MODULATED FUNCTION 
AI - IMAGINARY PART OF MODULATED FUNCTION 

MR - DUMMY_COMMON ARGUMENT _ _ 

MI “ DUMMY COMMON ARGUMENT 


SUBROUTINES REQUIRED 
MAGTAP 


************* ^^** » * * » » *_* ***** * * * * 

DIMENSION ARC2S6) »AIC 256) #MRC 2048)# MIC 2048) " 

COMMQN/STORGE/ AR # AI » MR iM l 

11 = 1 

12 = 2 _ _ 

14 = 4 
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CaLL MAGTAPCI4»NR»AR»IUNT1) 

CALL MAuTAPt I4»NR»AR» IUNT2) 

C*** CuNPUTE EXPONENTIAL FOR EACH COMPONENT 
LO 4 J=1»NC 

CALL MAGTAPL 11 » NR »mR»IUNT1) 

CALI MAGTAPE II »NR» Al » IUNtl) 

LO 3 1=1 * NR ' 

X=SETA*aRLI) 

ARC I ) =COSC X) 

AILI)=$1nEX) 

3 CONTINUE 

CALL MAETAPLI2»nR*AR»IUNT2) 

CALL MaGTAPC IE»NR# AI * IUNT2) 

4 CONTINUE 

CALL MAGTAPE I4»uR#AR» IU wTl) 

CALL MAGTAPE 14# NR »AR» IUn12) 
keTokn _ 

END 
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SUbKOUTlNE TSGEhC IUnT * MR* MC»DT* I OPT * IFM» A*F»THETA*KT »TP) 

CD* ********************************** 

C‘, 


C O 
CD 
CD 
CD 
CO 
CD 
CO 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CD 
CO 
CD 
CD 


PROGRAMMER and date 

RICHARD C. THOMAS 
TRW SYSTEMS 
DONE' 1972 


PURPOSE 

GtNERAIES THE bUILT IN TEST SIGNALS FOR FREQUENCY MODULA- 
TION! FM) OR PHASE MODULATION. THE TEST SIGNALS ARE SUM OF 
SINUSOIDS ’PER IODIC FOUR LEVEL TEST SIGNAL* AND A SQUARE 
WAVE. 

USAGE 

CALL TSGEuC IUNT » MR » MC » DT » IOPT r IFM # A F * THETA * K » TP » BET A ) 
DESCRIPTION OF PARAMETERS 
Input 

CALLING SEQUENCE 

IUNT - TAPE UNIT FOR GENERATED OUTPUT SIGNAL 
MR - NUMBER OF ROWS IN OUTPUT ARRAYS 

MC - NUMBER OF COLUMNS IN OUTPUT ARRAYS 

DT - DELTA TIME BETWEEN OUTPUT POINTS 

I OPT — *£Q. l SINUSOIDAL TEST SIGNAL 

- .EG. 2 THREE LEVEL GRAY TEST SIGNAL 

- .EQ.3 PERIODIC SQUARE WAVE TEST SIGNAL 

IFM - .EQ.l FM MODULATION 

•r .EQ.2 PHASE MODULATION 

a t .ioPT.Eaa amplitude 

F - . I OPT. EQ.l FREQUENCY 

THETA - .IOPT.EQ.l PHASE 

K - .IOPT.EQ.l NUMBER OF POINTS 

TP - . IOPT .EQ.2 DEFINES TIME POINTS OF BREAKS 
COMMON 

AR - REAL PART OF SIGNAL BEING GENERATED 
Ai - IMAGINARY PART OF SIGNAL BEING GENERATED 
WR - DUMMY COMMON ARGUMENT 

WI - DUMMY COMMON ARGUMENT 
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CL 

CL OUTPUT 

CD CAuL SEQUENCE 

Cu COMMON 

CL PRINT 

CD TAPE 

CD 

CL SUdKOU FINES REQUIRED 

CL MAGTAP 

CD 

CL* ********************************** 

DIMENSION AL2b) »FC25) »THETAt25) »TPC3) »AR1256) *AI1256) 

DIMENSION wRL 2048) »WlC2048) 

COMMON/ST LRGc/ AR » A I » wR » W I 
nR = MR 
NO “ i*iC 
(•I •• K 

TWOPiSb.263luS.31 
12 = 2 
14 = 4 

CALL MA0TAPCI4#NR»AR»IUNT) _ . 

C**** SET IFlAG 

JFLAG=I0PT+3 

C**** 1. TIALlZE 

COn-O.O _ 

DO 1 1=1 » NR 
. 1 Alt i)=u«0 

IFC IOPTtNE.l) GO TO 3 
DO 2 I=1?M 
FtI>=TwOPl*FCI> 

THET At 1 ) =0 . 0174532925*THETAC I ) 

2 CONTINUE 

3 0=1 . 

K=1 

T = U.U - ... 

S=0.0 

At 5) =0.0 .... 

ALAST=0.0 

tsum=o.u . .. _ 

smean-o . 0 
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5 CONTINUE 

1FL IFLAG.Eu.b) S=AL1) 

IFCIFLAG.NE.5) GO TO 90 

70 UO 71 1=1 » 4 

71 SMlAN=Si*.EAN+al I ) *C IPC5)-TPCI> ) 

SMEAN = SMEAN/TPCS) 

90 CQimT iNUE 

C**+* GO TO IFLAG OPTION 
10 CONTINUE 

GO TO L400»500#t>00) »IOPT 
C**+* SINUSOIDAL TEST SIGNAL - PHASE MODULATION 

400 DO 402 1=1 » NR 
T=L1-1)*0T + TSUM 
TEMP=0 . 0 

DO <*01 L=1»M 

401 TEMP=TEMP ♦ ALL)*SINCFCL)*T + THETACL) ) 

ARC 1 )=T tMP 

402 CONTINUE 
GO TO 999 

C**»* THREE GRAY LEVEL TEST SIGNAL - PHASE MODULATION 

500 DO 503 1=1 » NR 
IFLT.LT.TPU)> GO TO 501 
S=S + AU> 

K=K 4 1 

501 IFCK.NE.6) GO TO 502 
K=1 

T=C.O 

S=0.0 

502 ARL I )= S - SMEAN 

T= T ♦ DT~" 

503 CONTINUE 

GO TO 999 

C***» PERIODIC SQUARE tfAVE TEST. SIGNAL - PHASE MODULATION 

600 DO 602 1=1 * NR 
IFtT.LT.TPCD) GO TO jbOl 
T= T -TPC1) 

S= -S 

601 ARC I )=S 

T=T+DT 

602 CONTINUE 
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C**»* *KiTt TAPE 

999 CALL MAGTAPLI2HxR»mR#IUNT) 
CALL MAGTAPLI2»NR»AI*IUNT> 
TSUM=T + 01 

J= wl + 1 

iFtJ.Lt.NC) GO TO 10 
CALL MAGTAPU4»NR»AR» IUNT) 

return 

END 
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CD* * 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 


SUBROUTINE TTKAuSl 1A* IB» IUNITA* IUNITB) 

************************ 

PROGRAMMER AnD DATl 

RICHARD C. THOMAS 
TRW SYSTEMS 
JUNE 1972 

PURPOSE 

USAGE 

CALL TTRANSL IArIB»IUNITA»IUNIT3#MBUFP) 
DESCRIPTION OF PARAMETERS 


* * 


* * * * * 


* * 


CD 



CD 

input 


CD 

CALLING SEQUENCE 


CD 

IA - NUMBER OF WORDS/RECORD ON 

INPUT RECORD 

CD 

lb - NUMBER OF REAL DATA INPUT 

RECORDS 

CD 

IUNITA - INPUT DATA TAPE UNIT 


CD 

IUNITB - OUTPUT DATA TAPE 


CD 

MBUFF - SIZE OF WORKING BUFFER 


CD 

COMMON 


CD 

BJF - WORK AREA 


CD 

BUFI - WORK AREA 


CD 

BUFF - WORK AREA 


CD . 

BUFFI - WORK AREA 


CD 

TAPE 


CD 



CD 

OUTPUT 

■ 

CD 

CALL SEQUENCE . 


CD 

COMMON 


CD 

TAPE 

- . 

CD 



CD 

REMARKS AND RESTRICTIONS 



CD 

SUBROUTINES REQUIRED 


CU 

MAGTAP _ 


CD 




CD* ***********************'***** 

DIMENSION bUFL 256 KBUFFC 2048) »8UFFIC2048) »BUFK 256) 
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COmMuN/STO«oE/BUF » uUF I » BUFF » BUFF I 

II = 1 
Id = £ 

14 = + 

mbuff = a04d 

ibr ~ I«*lb 

1PI IbF .Ll .M ddFF) MoUFF = IBF 

C *** tOMPUT MAXIMUM ImUMdLR OF OUTPUT RECORDS TO BE LOADED PER COMPLETE 
C *** PASS OF INPUT TAPE. 

IX s HBUFF/IB 

C *** COwPUTc NUMBER OF COMPLETE PASSES OF INPUT TAPE. 
hAxRD = IA/IX 

III = IA/IA 
XY = IA/MAaKd 
11 = Ib/MAXRU 

DO 200 1MX = 1 » LIZ 
DO 100 IXX = 1»MAXRD 
DO SO L = 1 • ll 

CALL MAoTAPCIi. IA.BUF. IUNITA) 

CALL MAGTAPC II » 1A »dUFI r IUNITA) 

DO bu 1-1 » I Y 

OK s ir*ClMX-l) 4 I 

LL = L ♦ IZ*CIX V -1) + IZ*MAXR0*CI-1> 

SUFFCLi.) = dUFL jK) 

BUFFICLL) = BUFICUK) 

50 CONTINUE 
.100 CONTINUE 

DO 76 1=1. AX 

MX = Id*C I**l ) ♦ 1 

CALL rtAGTAPCI2»iB»BUFFtMX)»lUNITB> 

CALL MAGTAPC IZ» IB’BUFFICMX) » I UNITS) 

75 CONTINUE 

CALL MAGTAPC 14 > 14 »dUF» IUNITA) 

200 CONTINUE 

CALL MAGTAPC 14 >I4»BUF»IUNITB) 

900 RETURN 

END 
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CL* 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 


SUuROUTiNE TrfkITtt lUt NCtNR) 

********************************* 

PhDGKAwfiER AND DATE 

RICHARD C. THOMAS 
TRW SYSTEMS 
DUNE 1972: 


PURPOSE 


TH15 SUBROUTINE WILL RtAb THE REAL AND IMAGINARY PARTS 
A SIGNAL ON TAPE AND WILL OU1 PUT THE SIGNAL ON THE 

printer. 


OF 


USaGl 


CALL T#RITECIU»NC*NR> 
DESCRIPTION Of PARAMETERS 
Input 

CALLING SEQUENCE 


CD 

1U 

- 

DATA TAPE UNIT 


CD 

NC 

- 

NUMBER OF COLUMNS IN OATA 

MATRIX 

CD 

CD 

NR 

COMMON 


NUMBER OF ROWS IN DATA MATRIX 

CD 

AK 

• 

INPUT/OUTPUT BUFFER 


CD 

AI 

• 

INPUT/OUTPUT BUFFER 


CD 

WK 

• 

DUMMY COMMON VARIABLE 


CD 

CD 

Wl 

OUTPUT 

* 

DUMMY COMMON VARIABLE 


CD 

CD 

PRINT 




CD 

CD 

CD 

SUBROUTINES 

NOwE 

REQUIRED 


CD* 

******* 

* * * 

************* 

♦ * * 


* * 


DIMENSION ARC 256) »AIC256) »WRC2Q48> ' WIC2048) 

COMMON/ ST OKGE/ AR»AI»WRmkI 

11=1 

14=4 

CALL MA6TAPC I4r I4»AR# 1U) 
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DO bo u=lMC 

CALL MAGTAPt I1 »nR»AR» IU) 

CALL MAGTAPC11*NR»a 1»IU) 
NKlTtLb’2005) J 
20Ub FORMaTL ih » tiHCOLUMw = »I4) 

*RiTtto»2010> L mRC I ) » 1=1 rNR) 

2010 FOKMATl 1H »HHR£AL/Lt>t;i6.8) ) 
*KITtLo»20Ab) LaICI) »I=1»NR) 

201b FORMAT' 1H t 9HIMAGINARY/1 8E16.8) ) 
bO CONTI Nut 

CALL MAOTAPL I4» I4» AR» IU) 

RETURN 

tNU 
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Cu* 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD* 


jDDKODT iNE TIMEnl i »«J»M*S) 

********************************* 


PROGRAMMER AND DATE 

' IChARD C. THOMAS 
TRW SYSTEMS 
JunE 1<*72 

PDKPDSt. 

USING TWO CALLS TO SUBROUTINE CLOCK AS INPUTS TIME CALCU- 
LATES THE RUN TIME IN MINUTES AND SECONDS OF A SEGMENT OF 
THt PROGRAM. 


USmGE 

dESCKIPI ION of parameters 
IhPUT 

calling sequence 

I - INITIAL TIME 

J - FINAL TIME 

M - MINUTES 

S - SECONDS 


output 

CALL SEQUENCE 
REMARKS AND RESTRICTIONS 
SUdKOUl INES REQUIRED 
NONE 

********************************* 
X = J - I 
T = X/1024. 

M = T/6U. 

Y = M 

S = T - Y*60. 

RETURN . . . . 

END 
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CD* 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 


SUBROUTINE ISARLIUNIT1»IUNIT2»MR2»NC2»DT1»NR1»NC1»KSR) 

********************************** 

PROGRAMMER AND uATe. 

RICHARD C. THOMAS 
TRW SYSTEMS 
JUNE 1972 

PURPOSE 

ACCEPTS The USER SUPPLIED TEST SIGNAL FROM TAPE. RESAM- 
PLES IT AT A HIGHER DATA RATE AND GENERATES AN OUTPUT 

signal, tape compatible with sap. 


CALL ISARCIUNIT1»IUNIT2»NK2»NC2»DT1»NR1.NC1»KSR»BETA) 

description of parameters 
Input 

CALLING SEQUENCE 

IUNIT1 - INPUT TAPE UNIT 
IUNIT2 - OUTPUT TAPE UNIT 

NR2 - NUMBER OF ROWS IN OUTPUT ARRAY 

NC2 - NUMBER OF COLUMNS IN OUTPUT ARRAY 

DU - SAMPLING INTERVAL FOR INPUT DATA 

NR1 - NUMBER OF ROWS In INPUT ARRAY 

NCI - NUMBER OF COLUMNS IN INPUT ARRAY 

KSR - DEFINES SAMPLING RATE FOR RESAMPLED SIGNAL 

BETA - MODULATION INDEX. MULTIPLIES ENTIRE TEST SIGNAL 

COMMON 

AR - WORK AREA 

AI - «ORK AREA 

S - INPUT BUFFER 

SX - INPUT BUFFER 

TAPE . . - 

OUTPUT . . 

CALL SEQUENCE 

COMMON 

CARD 
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TAPE 

SObKOUT IftES REuUlrtbG 
fcAoTAP 

NONt 

'**»*«*****************« 

HMENSION AKC25b) f A1C256) »SL204fl) »SXC204d) 
cOMhON/STORGE/ AK » A I » S » SX 
CALL MAol APt K4 # hRl » S r XUNlTl > 

CALL i‘,AGTAPtK4»NRl»S»IUNXT2) 

Kl = 1 
*2 = 2 
*4 = 4 

FLTKSR = KSR 

CALL MaGTAPCK1»i>iR1»S»XUnITX) 
l> 0 b I=x»NH2 
AICX) = 0.0 
& COWTInul 
PI = SCi) 

P2 = SC 2) 

XI = 2 
12 = 1 
Jl = 1 • 

J2 = 1 

TERM = CP2-P1) /FLTKSR 
K = 1 

20 CONTXiWL 

ARC 12) = PI t TERM*FL0ATCK“1) 

12 = X2 ♦ 1 
K = K + 1 

1FCX2.LL.NR2) oO TO 40 
CALL M ACT APC K2 » NR2 * AR » I UN I T2 ) 

CALL MAGTAPC K2»NR2» AI » IUNl T2) 

X2 = 1 
u2 = 02 ♦ 1 
XFCJ2.Lb.NC2) GO TO 40 
CALL MAGTAPCK4'hRl'AK'IUNlTl) 

CALL MAGTAPCK4'NR1'AR'IUNIT2) 


SAP PROGRAM LISTING 


KtTUKN 
4b CONTINUE 

1FLK.LE.KSK) GO TO 20 

N - 1 

11 = II + 1 

IFLll.GT.NKi) GO TO GO 
ob CONTINUE 
PI = P2 
P2 - SLID 

TERM = LP2-PD/FLTKSR 
GO TO 2b 
O0 CONTINUE 
J1 5 Jl ♦ 1 

IFCdl.6l.NCl) Go Tu 100 
CALL KA6TAPCK1 hmR1»S»IUNIT1) 

11 = 1 
GO TO GO 
10b CONTINUE 

wR1TeC6’20U0) 

2000 FORMAT! b3H * * * NOT ENOUGH INPUT DATA TO GENERATE OUTPUT * * *) 
STOP 
END 
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SUBROUTINE MAgTAPCD 

cl* ********************************** 


CL 




CL 

PROGRAMMER AND 

GATE 


CL 

RiCHARu C. THOMAS 


CL 

TKv» SYSTEMS 


CL 

LOnE 

1972 


CL 




CL 

PURPOSE 



CL 

MAGTAP IS A LUMMY ROUTINE THAT ALLOCATES CORE STORAGE FOR 

CL 

A CS“ 

1 ASSEMBLY LANGUAGE ROUTINE NAMED MAGTAP. 

THE CS-1 

CL 

PROGKAM PEFORMS ALL TAPE OPERATIONS REQUIRED BY 

THE PRO- 

CL 

GRAM* 

FOR TAPE PARITY ERRORS THE MAGTAP ROUTINE 

ATTEMPTS 

CL 

RECOVERY A MAXIMUM OF TEN TIMES 


CL 

FOR UNRECOVERABLE PARITY ERRORS* FRAME COUNT ERRORS AND 

CL 

TIMING ERRORS THE COMPUTER COMES TO A HALT WITH 

NO PRINT- 

CL 

EL ERROR MESSAGES. THE OPERATOR IS ABLE TO DETERMINE THE 

CL 

NATURE OF THE ERROR faY INDICATOR LIGHTS ON THE 

TAPE DRIVE 

Cl 




CL 

usage 



CL 

CALL 

MAGTmPC IOPT * NUM # ARRAY * IUNT ) 


CL 




CL 

DESCRIPTION OF 

PARAMETERS 


CL 




CL 

INPUT • 



CL 

CALLING SEQUENCE 


CL 

IOPT 

- »EQ* 1 READ TAPE RECORD 


CL 

I OPT 

- .EQ.2 WRITE TAPE RECORD 


CD 

IOPT 

• *EQ*3 BACKSPACE RECORD 


CD 

IOPT 

- .EQ.4 REWIND TAPE 


CL 

IOPT 

- «EQ»5 WRITE FILE MARK 


CL 

IOPT 

- »EG«6 SPACE FILE FORWARD 


CD 

IOPT 

- »EQ»7 SPACE FILE BACKWARO 


CL 

NUM 

- NUMBER OF WORDS IN OATA RECORD 


CD 

ARRAY 

- LATA RECORD 


CL 

IUNT 

- TAPE UNIT 


CD 

TAPE 

- 


CD 




CL 

OUTPUT 



CD 

TAPE 
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CD 

CD KtMARKS AND RESTRICTIONS 

CD WHEN Tape operation is neither a read or write? num and 

CD array are dummy parameters. 

CD this routine allocates the necessary core required for 

CD The CS-1 program which must be loaded by paper tape, see 

CD PROGRAM OPERATING INSTRUCTIONS FOR CORRECT PROCEDURES. 

CD SUdROUTInES RLwUIRED 

CD NONE 

CD 

CD* ********************************** 

DIMENSION At 300) 

At I) = 1.0 

RETURN 

END 
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SAP PROGRAM LISTING 


SUbKOUT lNL TRFNLFREQ* IOPT »NORDER»FCEN*FCUT »RIP*TMAG»F»NF»NS»TR#TI 


CD* 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

Cu 


A»TPH) 

********************************* 

PROGRAMMER AND DATt 

RICHARD C* THOMAS 
TRW SYSTEMS 
DUNE 1972 

PURPOSE 

CALCULATES THE REAL AND IMAGINARY PARTS OF THE DESIRED 
TRANSFER FUNCTION 


* 


CD 

CD USAGE 

CD CALL rKFNtFREQ#IOPT»NORDER»FCEN»FCUT»RlP»TMAG»F#NF»NS» 

CD TR»T1»TPH) 

CD 


CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 


DESCRIPTION OF PARAMETERS 
INPUT 

CALLING SEQUENCE 

FREQ - FREQ AT WHICH TRANSFER FCN IS TO BE EVALUATED 
IOPT - .EQ.l BUTTERWORTH FILTER 

- .EQ.2 CHEBYCHEV 

- .EQ.3 INPUT DATA FILTER 

NORD - ORDER OF THE BUTTERWORTH OR CHEBYCHEV FILTER 
FCEN - CENTER FREQUENCY FOR BUTTERWORTH OR CHEBYCHEV 

FCUT - CUTOFF FREQUENCY FOR BUTTERWORTH OR CHEBYCHEV 

RIP - RIPPLE FOR BUTTERWORTH FILTER 

TMAG - MAGNITUDE OF INPUT FILTER TRANSFER FUNCTION 
F - FKEQUENCYCIN HERTZ) AT WHICH MAGNITUDE AND PHASE 

NF - DIMENSION OF TMAG»TPH AND F 

TPH - CORRESPONDING pHASECIN RADIANS) 


CD* ************ » » » ************_******* 
DIMENSION Ft 1 ) * TMAGC 1 ) * TPHC i f 


DIMENSION CC6*6)*CCC6) ... 

DATA CLCCIrJ)* 1=1*6) *J=1*6)/1. 0*0. 0*0. 0*0. 0*0. 0*0.0# 
11. 414213662*1. 0*0. 0*0. 0*0. 0*0.0* 


22 . 0 * 2 . 0 * 1 . 0 * 0 . 0 * 0 . 0 * 0 . 0 * 
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SAP program listing 


52. 615125929# 3. 414213562»2. 613125929# i.0»0.0»0.0» 

43. 2360o7977# 5. 2360o?977» 5. 236067977 #3. 236067977 #1.0 #0.0 
53. 663703306# 7. 464101615# 9. 141620172# 7.464101615# 
o3. 665705305# 1 .0/ 

IFtFREG.LT. 0.0) GO TO 800 
GO TO C100#200»300) # IOPT 
C*** oUTTERwoRTH FILTER 
100 CONTINUE 
IRS*v=-l 
IISW= 1 
FSAVE=1.0 
1SW=-1 
CR-1.0 
CISO.O 

FNORM= AttSC FREg-FCEN ) /L FCUT-FCEN ) 

LO 10 I=l#NURbtK 

FPROO“-FSAVE*FNORM 

X= CCl»NORL)ER)*FPROD 

1FCISW.LT.0) GO TO 3 __ 

tR=CR ♦ IRSw*X 

iSws-1 

lRSfc^-IRSIw 

GO TO 10 

3 CI= Cl + IISW*X 
ISw=l 

I1SM=-11SM 
10 FSAVE = FPROD 
b=CR**2 ♦ CI**2 
TR= CR/b 
TI=-CI/0 
RETURN 

C*** CHEBYSHEV FILTER 
200 CONTINUE 

2FC NS.GT.l) GO TO 201 

X=l. 0/RIP _ __ 

SINHI=AL0GL X+SGRTC X**2+l ) ) 

XORUER-NOROER 

A=SlNHI/XORuER 

TEMP1=EXPCA) 

TEMP2=EXPC-A) 



SAP PROGRAM LISTING 


COSh=L TEMP1+TEMP2 ) /2. 0 
SlNh=L TEMP1-TEMP2 > /2 . 0 
TANH=SIwH/COSH 
201 CUNT iNUE 

FNORM= ABSL FREO-FCEN ) /C FCUT-FCEN ) 

GO TO L210»220»230»220r250»220i»NORDER 
210 CCLA)=1.0/TANH 
GO TO 269 
220 CONTINUE 

r,RiTtCb»2000) NGRDER 

2000 FOKMATL 1H1 * 20X90HYOU ARE IN A HEAP OF TROUBLE * BOY. NORDER HAS SEE 
IN ALTERED BY PROGRAM TO AN INVALID VALUE* 12) 

STOP 

230 CONTINUE 

TEMP = fANN**3 + 3.0*TANH 
CCC1) = C5.0*TANH**2 + 3.0) /TEMP 
CCL2) = 6. 0*TAWH/TEMP 
CCL 3) = 4.0 /TEmP 
GO TO 269 
250 CONTINUE 

TEMPI = TANH*0. 309017 

TEMP12 = TEmP1**2 

TEMP2 = TANH+0.b09017 

TEMP22 = TEMP2**2 

TEMP3 = TEMP12 + 0.9045085 

TEMP4 = TEMP22 4 0.3454915 

TE.MP5 - TEMP,5*TEMP4 

CCC1) = 1.0/TANH + 2 . 0*TEMP1/TEMP3 + 2.0*TEMP2/TEMP4 
CCC2) =C2.0*TEMP1/TANH)/TEMP3 ♦ C 2 • 0*TEMP2/T ANH ) /TEMP4 ♦ 

1LTEMP3 ♦ TEMP4 ♦ TANH**2) /TEMPS 
CCL3)=1.0/CTANH*TEMP3) ♦ 1.0/C TANH*TEMP4) ♦ CCTANH + 2.0*TEMP1 
1 ♦ 2.0*TEKP2>/TEMPS> 

CCC4) = C2*0*L TEMPI ♦ TEMP2) + TANH)/CTANH*TEMP5) 

CCL 5) = 1.0/C TANH*TEMP5) 

GO TO 269 _ . 

269 CONTINUE 
IRS*=-1 
1I5W= 1 

F SAVE=1.0 ... 

ISw=-l 


SAP PROGRAM LISTING 


CR=1.0 

1 1 = 0.0 

go £71 I=l»NORJER 
FPRCO=FSaVE*FNORM 
X=CCCl)*FPKOD 
IFLISw.LT. 0) GO TO 272 
CH=CR + IRSw*X 
1SW=“1 
1KSW=-1KSW 
GO TO 271 
272 CI=CI+I1SW*X 
IS»=1 

1ISW=-USW 
271 FSAV£ = FPROD 
0=CR**2 ♦ CI**2 
TK= CR/l 
T l=-Ci/U 
RETURN 

C*** INPUT DATA FILTER 

300 CONTINUE 

DO 301 N=NS»50 

IFCFREQ.lT.FCM-1)) GO TO 302 

301 CONTINUE 
NS=1 

GO TO 300 

302 NS=N 

DF=LFREu - FCNS))/CFCNS+1) - FCNS>) 
TMAGPsTNAGCNS) ♦ DF*C TMAGtNS+i ) - TMAGCNS) ) 
TPHP=TPHCNS) ♦ OF*LTPHCNS+l) - TPHCNSJ) 
TR=TMAGP*COSCTPHP> 

TI=TMAGP*SINCTPHP) _ 

RETURN 
GOO NS=1 
RETURN 
ENU 


3. PLOT GENERATION PROGRAM (PLTGEN) 


3.1 MODULE/SUBROUTINE DESCRIPTIONS 

The following list describes all the PLTGEN program modules and 
subroutines. 


ANN0TE 

Scaling and labeling routine 

C0MAND 

Constructs all plotter commands 

GRID 

Constructs plot grid 

INITAL 

Initializing routine 

MAGTAP 

Tape operations routine 

NUMBER 

Converts a floating point number to a Hollerith 
string 

PEN 

Pen select routine 

PL0TER 

Draws plot lines 

PL0TPT 

Repositions plot pen 

PLTGEN 

Main program for the plot software package 

PRINT 

Converts a string of Hollerith characters into 
the EAI plotter character >et. 


3.2 PLTGEN FLOW DIAGRAMS 

The flow diagrams for PLTGEN showing both the main program and all 
subroutines are shown in the following figures. 


3-1 



FLOW FOR MAIN PROGRAM PLTGEN 



SELECT PLOT PEN AND 
INITIALIZE SCALE FACTORS 


I 



CONSTRUCT GRID 



i 


1 SCALE AND ANNOTATE 1 

[THE GRID 




1 

r 

NORMALIZE AND SCALE 
AMPLITUDE DATA 

1 


PLOT SPECTRUM 

— 

i 


TERMINATE PLOT BY WRITING AN 

END OF FILE ON THE PLOT TAPE (UNIT 1) 

AND PRINTING "END OF PLOT" 

> 

, 


“T 
C » ) 


Figure 3-1. Flow for Main Program PLTGEN 
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FLOW FOR SUBROUTINE ANNOTE 



Figure 3-2. Flow for Subroutine ANNOTE 
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FLOW FOR SUBROUTINE GRID 



Figure 3-3. Flow for Subroutine GRID 


V'* 








FLOW FOR SUBROUTINE COMAND 



Figure 3-4. Flow for Subroutine COMAND 

(. 









FLOW FOR SUBROUTINE PLOTER 



Figure 3-5. Flow for Subroutine PLOTER 
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FLOW FOR SUBROUTINE NUMBER 



Figure 3-6. Flow for Subroutine NUMBER 
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FLOW FOR SUBROUTINE PRINT 


ENTER WITH 


X = X COORDINATE OF STARTING POSITION 

Y = Y COORDINATE OF STARTING POSITION 

ISTRNG = CHARACTER STRING ARRAY 1 HOLLERITH CHARACTER/ WORD 

NUMBER = NUMBER OF CHARACTERS IN ISTRNG TO BE PLOTTED 

IORIN = ORINTATION 0 = VERTICAL LETTERS, 1 = HORIZONTAL LETTERS 

DEL = DISTANCE BETWEEN CONSECUTIVE CHARACTERS 


EXAMINE A CHARACTER IN ISTRNG 



LOCATE EAI CHARACTER CODE 
AND CALL COMAND WITH SELECT 
CHARACTER COMMAND 


CALL COMAND WITH PLOT 
POSITION COMMAND 


INCREMENT CHARACTER 
SPACER BY DEL 


CHARACTER? 


RETURN 


Figure 3-7. Flow for Subroutine PRINT 







FLOW FOR SUBROUTINE PLOTPT 



Figure 3-8. Flow for Subroutine PLOTPT 
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FLOW FOR SUBROUTINE INITAL 



Figure 3-9. Flow for Subroutine INITAL 
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FLOW FOR SUBROUTINE PEN 



Figure 3-10. Flow for Subroutine PEN 
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3.3 PLTGEN LISTINGS 

i 

This section presents a complete listing of the PLTGEN program. 


O 


c 
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PLTbEN PROGkAM LISTING 



MAlNPKOoRAi'i PLToEN 


' 

CL* 


* * * 

**************** 

CL 




CL 

PROGRAMMER ANu LATt 



CL 

G. ... HOUSER 



CD 

TM SYSTEMS 



CL 

JUNE iy72 



CD 




Cl 

PURPOSE 



CL 

1HE PlIGeN PROGRAM 

i IS DESIGNED TO READ A DATA TAPE 

CL 

CONTAINING SPECTRUM DATA* 

GENERATED BY THE SAPOR I 

CL 

PROGRAM. TrlE SPECTRUM DATA IS SCALED AND CONVERTED INTO 

CL 

PLdT COMMANDS FOR 

THE LAI 

PLOTTER. A MAGNETIC TAPE 

CL 

CONTAINING ThE FORMATTED 

PLOT COMMANDS IS GENERATED. 

CL 




CL 

DESCRIPTION OF parameters 



CL 




CL 

luPul 



CL 

TAPE 



CL 

UNIT 2 C Ii.PuT LATA 

, TAPE) 


CL 

RECORD 1 C 36 WORDS) 

. 

CL 

WoRO 1 

IPLPOS 

- PLOT POSITION FLAG 

CL 

2-31 

ITITLE 

- TITLE Cl HOLLERITH CHAR. /WORD 

CL 

32 

DELTAF 

- DELTA FREQUENCY 

CL 

33 

FCENTR 

- CARRIER FREGUENCY 

CL 

34 

FSTART 

- STARTING FREQUENCY 

CL 

3b 

NPTS 

- NUMBER OF SPECTRUM POINTS 

CL 

3b 

IFC 

- FLAG FOR CENTERING CARRIER 

CL 

RECORD 2 L 480 

WORDS) 


CL 

WORD 1 - 480 

LATA 

- SPECTRUM VALUES CO.Q TO -70.0) 

CL 




CL 

OUTPUT 



CL 

CvMMOn 



CL 

labeled Common /a/ 

DATA C 480 ) - SPECTRUM VALUES 

CL 

CARD 



CL 

NONE 


. 

CL 

print 



CL 

prints »ENU of plot* when 

1 ALL CALCULATIONS ARE COMPLETE 

CL 

TaPE 



CL 

Ui»1T 1 LOUTPUT PLOT TAPE) 



c 


PLTGEN PROGRAM LISTING 


CL 


RECORD 1 CALL 

RECORDS CONTAIN 16 WORDS) 

CD 


wOKD 1 



Cu 


BYTE l 

an 

X COMMAND code 

CD 


Z 

an 

THOUSANDS BICIT 

CD 


3 

- 

HUNDREDS DIG IT 

CD 


4 

- 

TENS DIGIT 

CD 


5 

m 

UNITS DIGIT 

CD 


WORD Z 


BLANK 

CD 


WORD 3 



CD 


BYTE 1 

- 

SIGN 

CD 


2 

- 

NOT USED 

CD 


3 

• 

BLANK 

CD 


4 

- 

blank 

CD 


5 

- 

BLANK 

CD 


WORD 4 - 

12 

BLANK 

CD 


WORD 13 



CD 


BYTE 1 

• 

Y command code 

CD 


2 

mm 

THOUSANDS DIGIT 

CD 


3 

- 

HUNDREDS DIGIT 

CD 


4 

- 

TENS DIGIT 

CD 


.5 

- 

UNITS DIGIT 

CD 


WORD 14 


BLANK 

CD 


WORD 15 



CD 


BYTE 1 

• 

SIGN 

CD 


2 

• 

NOT USED 

CD 


3 

- 

BLANK 

CD • 


4 

• 

BLANK 

CD 


5 

- 

BLANK 

CD 


WORD 16 


BLANK 

CD 




• 

CD 

REMARKS 

AND RESTRICTIONS 


CD 

NONE 




CD 





CD 

SUBROUTINES REQUIRED 



CD 


ANNOTE 


- 

CD 


GRID 



CD 


inital 



CD 


MAGTAP 



CD 


PEN 



CD 


PLOTER 




C 
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o n n o o non non 


pltgen program listing 


cu plotp r 

Cu 

CO* ******************** -v ************ 

DIMENSION RdCORL/C 36) » NW0RDC2)» 1TITLEC30) 

COMMON /A/ OaTAC^oO) 

DATA ldLANK/5H / 

10=0 

11=1 

12=2 

14=4 

I5=b 

16=6 

17=7 

A36=36 

X9=-9.999 


RtrtINu TAPES - PLOT TAPE = UNIT 1# OATA TAPE = UNIT 2 

CALL MAbTAPtI4f 12»oATA»I2) 

KEAD 1ST DATA RECORDC 36 WORDS) FROM UNIT 2 

CALL MAGT API 1 1 r A 36 r RECORD » 1 2 ) 

ISX=2222 

1ST=3333 

ibx=0 

a6Y=0 

1DX=0 

1UT=0 


IF IPuPOS = 0 UPPER PLOT 
1 LOWER PLOT 


lPLPOS=RECOROtl) 
IFUPLPOS.Evl.l) IDY=9999 

UO b 1=1 » 30 
XTITLEL1)=HECOroLA+1) 
b CONTlNOc. 


non ooo nor onr on r r:nr. oor 


1 


I 

4 

$ 

1 

3 


PLTGE'N PROGRAM LISTING 


UELTAF=KECwftOL. 32 ) 

FCLNTK=KLCukl)L 33 ) 

FSTART=KEC0KuC34) 

uPTS=RlCQRuL3S) 

lFC=REC 0 RDL 3 b> 

initialize, plot routines 

CALL PELL II) 

CALL INlTAl.L ISa» 1 ST » IDXr IDY » IdX» 1 DY ) 
CoNSTkuCT toft id 

call GRID 

SCALE * LAt>EL» hND ANNOTATE ThE GRID 
CALL AnNOTEC IT! TLE»FSTART »DELTAF »FCENTR* IFC ) 
READ SPEC fftUM DATA FROM UuIT 2 


14ttU=4oC 

CmLL MAGTAPL II » 1430 » DATA » 12) 

NORMALIZE AND SCALE DATA FOR PLOTTING 
DO 1 U i=l»NPTS 

DATAC I ) =C L DATAC i ) ♦ 70 . ) /70 . ) *7 . + 1 . 

XO CONTINUE. 


PLOT SPECTRUM 
CALL PlOTERCnPTS) 

WRAP-UP PLOT BY. POSITIONING PEN AND REWINDING. TAPES 

CALL PEN Cl 7 i . 

CALL PLQTPTLX 9 »X 9 > 

CALL MAGTAPLI 5 »U»OATA»Il> 

CALL MaGTAPE 14 » 12 *uATA» 12 ) 


'k 

T 



PLT&EN PROGRAM LISTING 


HRITEL b r 2) 

2 FORMATL 12H END OF PLOT) 
STOP 


PLTGEN PROGRAM LISTING 


SUBROUTINE AiMNOIEL IUTLE*FSTART*DELTAF*FCEnTR» ICF) 

CD* ********************************** 

CO 

CO PROGRAMMER AND DATE 
CO G. L» hOUSEK 

TRW systems 

JO.mE 1972 


SUBROUTINE ANNOTE IS USED TO SCALE AND LABEL THE X AXIS* 
LABEL THE Y AXIS* AND PUT THE DESIRED INPUT VARIABLES l.M 
THE TITLE BOX. 


CALL ANNOTt C ITITl£*FSTART»DELTAF*FCENTR* ICF) 


CO 

CD 

LO 

Cu 

Cu 

CO 

CO 

CO 

CO 

CD 


PURPOSE 


OSAGE 


CO 

CO 

CD 

CD 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CD 

CD 

CD 

CD 

CD 

CD 

CO 

CO 

CD 

CO 

CO 

CO 

CO 


DESCRIPTION OK PARmMETERS 

input 

calling sequence 

ITITLE - TITLE OF PLOT* TO BE PLACED ON THE TOP LINE IN 
THE TITLE BOX. 30 CHARACTERS MAXIMUM. 

FSTART - STARTING FREQUENCY IN HZ. 

" ' OELTAF -‘DELTA FREQUENCY* INCREMENT BETWEEN CONSECUTIVE 

FREQUENCY LINES IN HZ. 

FCENTR - CENTER FREQUENCY OR CARRIER FREQUENCY IN HZ. 

ICF - FLAG FOR DETERMINING IF CENTER FREQUENCY IS 

CENTERED IN CENTER OF GRID. CIFC = 0 - CENTERED) 

output 

NONE 

REMARKS AND RESTRICTIONS 

ITITLE MUST BE AN ARRAY HAVING ONE HOLLERITH CHARACTER 
PER tftORD. THEflE 1$. A MAXIMUM OF 30 CHARACTERS IN ITITLE. 

SUBROUTINES REQUIRED 
PEN 

PuOTPT 

print 



METhUD 


PLTGEN PROGRAM LISTING 


NUMBER 


Hit LEnTEK FREQUENCY bias is subtracted from the starting 
FREQUENCY TO OBTAIN THE ACTUAL PLOTTED STARTING FREQ. 

Tht MAXIMUM FREQ. IS EXAMINED AND A SCALE FACTOR IS 
DtTERMlNEw. ThE POSITION AND VALUE OF THE FIRST TIC MARK 
IS CALCULATED AND THE DISTANCE BETWEEN TIC MARKS IS 

Found, the tic marks are placed on the x axis, the tic 

MARKS ALONG THE X AND Y AXIS ARE THEN LABELED. THE TITLE 
AND OTHER DESIRED PARAMETERS ARE PRINTED IN THE TITLE 
BOX THE STARTING FREQ. IS OMITTED FROM THE TITLE BOX 
ViHEN 1FC = 0. 


************** 

DIMENSION LABXLdO) »LAtJYC3Q) 
DIMENSION IhERTZt-4) 
DIMENSION ICFU30) 


******************* 
f ITITLEL 30 ) ♦ ITLEC 30 > » IFSTt 30 ) » IDELFC 30 ) 


DATAC IFSTt I ) » I=1»2U) /5HS 

» SHT 

»5HA 

»5hR 

»5HT 

9 

Shi »SHN »5HG 

»5H 

»5HF 

»5HR 

»5HE 

9 

5HQ * SHO t SHE 

5H- / 

f 5HN 

»5HC 

»5HY 

»5H 

9 

UA 1 At ILF Tl I ) » 1=1 » 20 ) /5HC 

1 5HA 

»5HR 

»5HR 

»5HI 

9 

5Hl • * SH.R t SH 

»5HF 

» SHR 

»5HE 

»5HQ 

9 

bHD t SHE t 5HN 

5H= / 

»5HC 

»SHY 

*5h 

»5H 

9 


DATAL IDELFC I ) » 1=1 i 20 ) /5H0 

» 5HE 

»5HL 

»5HT 

» 5HA 

♦ 

Sri » ShF » 5HR 

» she 

r5HQ 

» SHU 

» 5HE 

♦ 

SHN »BHC »5HY 

1 5H 

»5H 

1 5H 

1 5H 

♦ 

5H= / 





DATA 

CLABYLi) »I=l»iS)/5hP 

» 5HO 

» 5HW 

»5HE 

»5HR 

♦ 

5H» 

»5H 

1 5HD 

»5HE 

»5HC 

* 

SHI 

*5HB 

»5HE 

* 5HL 

* 5HS 

DATA 

CLABXti) rI=l»ll)/5HF 

* 5HR 

»5H£ 

» 5HQ 

» 5HU 

* 

5HE 

»5HN 

»5HC 

»5HY 

»5H» 

♦ 

5H 

/ 




DATA 

ClhEKTZCl) »I=1»4)/5HK 

» SHM 

r5HH 

»SH2 

/ 


DATA IBLANK/5H 
INTEGER DP r DOWN 
DELF=OtLTAF 
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non o on n n r r> 


plTgen program LISTING 


FSTk=FSTART 

FCNTK=FlENTR 

FCNT=FCImTR 

LlL rtRMlNL SCALED CBIASED) STARTING FREQUENCY dl 

SUBTRACTING OFF CARRIER FREQUENCY 

FSTKT=FSTART - FCNTK 
FCnTr-0 • 0 
IFC=ICF 
UP=7 
00*N=6 
uELX=*Uobb 
10=0 
11=1 
12=2 
114=14 
115=15 
120=20 
130=30 


CALCULATE STARTING FREQUENCY IF IFC=0 

IFL IFC . LQ . 0 ) FSTRT=FCNTR-240 . *DELF 

CHOOSE SCALE FACTOR FOR FREQUENCY LABELS 

FMAX=4b0 . *DElF+FSTRT 
FMAX=AMAXlt ABSCFMAX) » ABScFSTRT) > 

IFtFMAX-l.E+4) 20*2Q» 10 
10 IFLFMAX-l.E*?) 30»30»40 
"0 CONTINUE 

LABXL 12)=lHERT*L3> 

LABXL 13) =IH£RT2L 4) 

LAbXL 14)=IttLANK 
SCALE=l.E+00 
GO TO 50 
30 CONTINUE 

LABXL 12 )=IH£RTZC1) 

LAbXL 13 ) = IHERT2L 3 ) 





O 


PLTGEN PROGRAM LISTING 


i 

< 



LAbXL 14 i = I ;iER T ZL 4 ) 
bCALE=i »E+b 
GO TO SO 
40 COnTINUL 

LAflXLl2)=lHERTZl2) 

LAoXL lb) = lHtRTZL h) 
uAbXCl4)=lHLRTZL4) 

SCALt-1 .t+o 
bO CONTINUL 
C 

c determine fklguenct range and step size 

c 

fstrts=fstrt/scale 

DLlFS=UELF/SCAll 
FRANgE=hBG • *l»ElF s 
STEP=FKANGl/2*. 

1FLSTEP-1 * 0 ) t>u » 60 • 70 
oO CONT InUc 
ST tP=.b 
GO TO 100 
70 CONTINUE 

IFLSTEP-10.) 74 » lOu » 77 

74 CONTINUE 
ISTEP=STEP 
bSTEP=ISTEP 
ISTEP=ISTEP+I 
IFCSSTEP-STEP) 7b»75»76 

75 ISTEP=ISTEP-1 
7 b STEP=1STEP 

GO TO 100 
77 CONTINUE 
1=0 

60 CONTINUE 
1=1 + 1 

IREto=STtP/lO.**I 
IFL1REM.LT. 10) GO TO 90 
GO TO 60 
90 CONTINUE 

STEP=C IREM+1 ) *10**1 
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n n n n n n 


PLTGEN PROGRAM LISTING 


DETERMINE POSITION OF FIRST TIC MARK AND FIRST FREQUENCY 

number 

100 IF lKiT=FSThTS/STEP 
STARTY=1. 

STARTX=-6.222 

XSluP=7.777 

DX=.03j33333/uELFS 

FlRST=iFIRSr*STLP+STEP 

F1KSTX=STARTX+' c IRSr-FSTRTS)*DX 

DELTAX=uX*ST£P 

PUT TIC MARKS ON X GRID 


TIC=.0t> 

TICY=STARTY-TIC 
TIC2=STARTY-2.*TIC 
CALL PEUlUP) 

CALL PLOTPTC STARTX t ST ARTY ) 
CALL PENCDOMN) 

CALL PLOTPTCSTAKTX» TIC2) 
CALL PEnCUP) 

FIRS T2=F IRSTX-DELTAX/2 . 
X=F1RSTX 
X2=F IRST2 
110 CONTINUE 

IFCX-XSTOP) 120 1 120» 130 
120 CONTINUE 

CALL PL01PTCX2»STAkTY) 

CALL PENCDOwN) 

CALL PLOTPTC X2#l ICY ) 

CALL PENLUP) 

CALL PLOTPTC X»STARTY) 

CALL PENLDOnN) 

CALL PL0TPTCX»TIC2) 

CALL PENC UP) 

X=X+DELTAX 
X2=X2+D£lTAX 
GO TO 110 
130 CONTINUE 



non non r> n n 


PLTGEN PROGRAM LISTING 


LABEL TIC MARKS WITH SCALEO FREQUENCIES 

CAuL PENCIL) 

Y-.76 

A=STARTX-.co 

CALL NUMBER! FSTRTS* 10 »X»Y) 

FRACT=.000l 

X = FIRSTX~.<£ 

COUNT=FlRST 
140 CONTlNUt 

IFtX-XSTOP) 14S» 145» 170 
14b CONliuUfc 

XFISTEP-1.0) 160 » 160 » 150 
ISO ICOUNT=COUNT+FRACT 

1FL ICOUuT.LT . 0 ) ICGUNT=COUNT-FRACT 
C0UNT=1 CUUnT 
lbO CONTINUE 

call NUMBER! COUNT » 1 0 # X » Y ) 

X=X+u£lTAX 

count=count+step 

GO TO 140 


PRINT X LABEL AT CENTER BOTTOM OF PLOT 


170 CONTINUE 

00 180 i=l»30 
ITLE! I ) = IT I TLEC 1 ) 

180 CONTINUE 
DY=1. 
x=-.s 
Y=.l 

CALL PRINT! X » Y » LABX r 1 14 1 1 1 r OELX ) 

X=-9.200 

Y=7.5 

UELY=.i 

CALL PRINT! X» Y»LABY » I15» IO*_DELY> 

LABEL Y AXIS WITH DECIBEL SCALE VALUES CO.O 
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Y=tt* J 
X=**9.Qt> 

AN=1U. 

00 190 1=1 * 8 
XN=XN-1G* 

CALL NUMBER! XN» 10 *X*Y) 
Y=Y-UY 

190 CONTINUE 


PRINT TITLE* CENTER FREQ* r AND STARTING FREQ IN TITLE t>0X 


XT=5*bGo 

YT=9.333 

CAlL PRlN’t l XT# f T * ITLE# 130* II *DELX) 
YT=YT-*3 

CALL PRINT! XT * YT * IOELF *120*11# DELX > 

X=XT+21.*0LLX 

CALL NUMBERC ULLF » II * X * YT) 

YT=YT-.3 

CALL PR1nT1XT»YT»ICFT#I20»I1#DELX) 
CALL NUMBERCFCNT*11*X»YT) 

OMIT STARTING FREQ IF IFC=0 

IFClFC.EQiU) GO TO 200 
YT=YT-*3 

- CALL PRINT! XT » YT » IFST # 120# II* DELX) 
CALL NUMBER! FSTR* Il*X» YT) 

200 CONTINUE 
RETURN 
END 
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SCuKoLliNL GRiu 

UL-* ********************************** 


Cw 

wl/ 

CL 

iu 

CL 

Cu 

CJ 

•CL) . 

CL 

Cl* 

CL 

LU 

CL 

Cu 

CL 

CL 

CL 

CL 

CD 

Cu 

CL 

CL 

CL 

CL 

Cl 

CL 

CD 

Cl 

CL 

CL 

CL 

CL 

Cl 

Cl 

CL* * 


rRugnA •*,. ,t.rt mwu wAlt 

G • u« nOLouK 
TR* Sf STEMS 
lone ‘ 1972 

PLkPujl. 

ThIs oLURlUTINl IS USED Tl DRAW A RECTANGULAR BOX 25 In. 
Lif 7 IN. '* 1 1 H TIC MARKS 1 INCH APART ALONG THE Y AXES. 

A SMALL TITll lOX IS ALSO DRAWN l.N THE UPPER RIGHT 
CLKWLR. lhuSL DIMENSIONS ARE VALID ONLY WHEN USING A 
l.o TO 1. A SCALE FACTOR ANU A 1. TO 1. Y SCALE FACTuh. 

LbrtC-L 

Ca„i_ OkIL 

LtSCKlPllOu OF PARAMETERS 
InPLI 

CALLlNo SEQUENCE 
noh£ 

KEMAKKb ANU RlSIRICTIONS 
None 

SooRlLTINES RluoIKEL 
PEN 

PLOTPT 

METHOD 

The KtCT ANGULAR box is drawn counter clockwise starting 
in THE LOnLR LEFT CORNER. ThE TIC MARKS AND TITLE BOX 
ARE DRAWN AS THEY ARE ENCOUNTERED. 

********************************* 

infEotR uP»DOW.'i 

iO=0 

11=1 



r n r nr 
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lt--c 

L,Wrth — O 

o»'-7 
I IC=.l 

bl HKTX=-6.bbb 
blAKTY=l . 

CALL PlNLUP) 

call PLOTPUsrAhUfbiHKTr) 
lmlL PlUCDOmn) 
a = S I hKTA+ • o33b 
l/LlI m— 1 • 

L 

L UkAw LOrtLK X AXIS LINE 

C 

uO lU 

LAlL PLOT Pi L X * ST AK f Y ) 
lKi.tu.ib) GO TO 10 
lKi.Lw.17) L»tLTA=.h333 

A=A+UtLl A 
lO LUNl I uUl 

C 

UK A a RIGHT Y AXIS wlTH TIC MARKS 

TlCX-X+1 IC 
l=SlAr\lt 
UV 20 l=i» 7 
l=Y+i. 

CALL PlOTPTLaK) 

IFLI.lO.7) GO TO 20 
CALL PlOTPTLT1CX»Y) 

CALL PLu‘ r PH.XK) 

20 COn I iNoc. 


DRAm TITLt BOX IN UPPER RIGHT CORNER 

YbOX=YK.b 

CALL PloTPKX» fbOX) 

XbOX=X 

CO Jo 1-1 tk 

auoX=XdwX“1.^ 



r r. r 
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CAuu PcOTPl L XdOX» YoOX) 
jU c0i«TIN0fc 

Call PLOT P 1 L AtJOA t Y ) 

CALL PLUtUr*) 

CALL PcGTPTlX»Y) 

CALL PLNlOOwN) 

UKAft UPPER X AXIS 

liCY=Y+TlC 
A-A” • 

uEcT A— 1 • 

L'O HU i = l»ld 
CAcL PcuTPTLXrY) 
if * l I « lu • IB ) GO TO HO 
1FL 1 «Ev»» 17) DEU A=.S333 
X-X-ULLIa 
40 COhUNUE 


C DRAW LEFT Y AXIS WITH TIC MARKS 

C 

TICX=X-i 1 C 
LO SO 1=1 » 7 
1 =Y- 1 . 

CALL PLOTPUXrY) 

IFLl.EO. 7 ) o 0 TO 50 
. call plotptlticx.y) 

CALL PlGTPTCa»Y) 

SO COnTINUE 

c 

RETURN 

END 
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Cu* 

cu 

CD 

Cu 

CU 

Cu 

cu 

cu 

cu 

cu 

cu 

CU 

cu 

Cu 

Cu 

Cu 

cu 

cu 

cu 

cu 

Cu 

Cu 

cu 

Cu 

cu 

cu 

cu ’ 

cu 

cu 

cu 

cu 

cu 

cu 

cu 

cu 

cu 

cu 

cu 

CD 


SUbhUUlilMt COMImU L iXCOM» IYCOM» 1/Xt IYY) 
********************************* 

PROGRAMMER ANU uATu 

G. L. HOUSER 
TRw SYSTEMS 
JU.mE 1972 

PURPOSE 

TO CONSTRUCT AND WRITE ON MAGNETIC TAPE# ALL EAI PLOTTER 
COMMAihUS* 

uSAGc. 

call COMAhU t IXCO m* IYCOM# IXX» IYY) 

UtSCKiP. IO h of parameters 
Input 

calling sequence 

lACOM - X command code 

1YCOM - Y COMMAND CODE 

IXX - X NUMERICAL VALUE TO BE PLOTTED. 

IYY - Y NUMERICAL VALUE TO BE PLOTTED. 

OUTPUT 

TAPE 

MAGNETIC TAPE UNIT 1 CSEE REMARKS BELOW) 

REMARKS ANU RlgTRICTION 

The On I VAC 1230 COMPUTER word CONSISTS OF 5. SIX 29T 
BYTES. EACH DECIMAL DIGIT C CHARACTER) IS REPRESENTED tjY 
ONE BYTE. THE FORTRAN MA6TAP ROUTINE MUST WRITE TEN 
BYTES FOR EACH 5 BYTE FORTRAN WORD SPECIFIED. THE FIRST 
S dYTES CONTAIN THE COMPUTER WORD SPECIFIED AND THE 
SECOND 5 BYTES ARE ALWAYS BLANK. 

THE EAI PLOTTER CAN BE MANUALLY SET TO READ RECORDS Anu 
WOKu SIZES OF VARIOUS LENGTHS. EACH PLOT COMMAND MUST 
HOWEVER CunTAIN 6 BYTES OF INFORMATION. THE PLOT CONSOL 
IS THEREFORE SET TO READ A 6 WORD RECORD# CONSISTING OF 



4 





CU 
CL 
CL 
CL 
CL 
CL 
CL 
CL 
CL 
C L 
Cl 
CL 
CL 
CL 
Cu 
CL 
CL 
CL 
CL 
CL 
Cu 
CL 
CL 
cu 
CL 
CL 
CL 
CU 
CU 
CD* 


PLTGEN PROGkAM listing 

ic BUc. IchaRACT c R) WORDS. THE FIRST 12 CHARACTER WORL 
Ito iHt X COMMAND word anl the sixth word is the y COMMAND 

a URL* 

bo in the a and r command *ords on the plot tape have the 

Same STRUCTURE WHICH IS AS FOLLOWSO THE X COMMAND IS THE 
FIRST bYTc. THE SECOND THROUGH THE FIFTH BYTES ARE THE 

thousandths* hundredths* tenths* and units positions of 
i H e value to be plotted, the sixth through the tenth 

btltS ARE NOT USEO. THE ELEVENTH BYTE IS THE SIGN ANU 
The TaELVTH aYTE IS NOT USEu. 



The InORu ARRAY IS CONSTRUCTED IN THIS ROUTINE. IT IS M 
SlKUCTUREu 10 SATISFY THE ABOVE REQUIREMENTS* KEEPING Xu m 

MIND That ihe kagtap routine writes a blank word after a 
each a character woru specified to be written on the mag.S 
tape. m 


SUohoUl INEb KLJUlKC.0 
MauTaP 


me I HOD 

Thc X COMMAuD code and numerical value are determined a 
The x PL01 COMMAND IS FORMED IN IwORDCl) AND IW0RUC2). 

'Iaokul7> anl irtORDtb) contain the properly formatted 
Y command, the 6 word record is then written on tape 
Unit I. tape unit 1 then actually contains A 16 word 
hECORU* LuE to the intelaced blank words. CSEE REMARKS) 

********************************* 

DIMENSION IisOKjL 6) t ICHARC 5) 

11=1 
12=2 
16= to 

IXSluN=c 
1Ys1GN=0 
1X=*XX 
IY=1 YY 
1XC=IXC0M 
1 YC=lYCuM 
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C 

0 Del EKi-ii Nti SIGN OF X AND Y NUMERICAL VALUES AND SET SIGN 

C *OKDS 

L 

IFllX.ei.O) IxSlGi!=HO 
IFCIY.LT.O) IYSIGN=4U 
AX=lMbSL 1 A 1 
AY=iAbsLXY) 

DO 10 I~1 1 AD 
iWoKDC I ) = 0 
1U CONTINUE 

c 

C LOAD X COMMAND CODE 

C 

IChAkC A)=1aC 
C 

c isolate each digit ok x numerical value 

c 

DO 20 A = lr4 
K = t>“ 1 

ITtwP=lX/lU 
IC;:AKf.N) = lX-ITEMP*10 
AA=1TLMP 
20 COnTINUL 

c 

C bUlLD X COMMAND fcORDS 

C 

IbORt <L 1 ) = ICHARL 1 ) *2**24+ICHARC2) *2**18+1CHARC 3) *2**12+ 

♦ ICHARC4) *2**6+ICHARC 5) 

1<*0RDC2)-IXSI6N*2**24 

C 

C LOAD Y COMMAND CODE 

C 

iCHAKll)=I*C 

C 

C ISOLATE EACH DIGIT OF Y NUMERICAL VALUE 

C 

DO 3b i = l »4 
K=6-I 

IILMP=IY/10 
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ICMAKLi\)=It-i rL.».P*lU 
IY=ITEMP 
3u CONTINUE 
C 

C BUILD V COMMAND WORDS 

C 

U.UkuC7> =lCHAKt i ) ♦ <i**24+ICHAKC2) *2**le+ICHARL3) *2**12+ 

* IChaHL4)*E**6+ICHARE5> 

iwORDta)=lYSiOU*2**24 
C 

C WHITE tIGHT FORTRAN WORDS ON UNIT 1 ESEE REMARKS ABOVE) 

C 

CALL MaGTAPL 12 i ib» IwOKDt ID 
C 

HE rt»HN 
tNU 
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SDdRoUIINE PuOTEK LNPTS) 

CD* ********************************** 


CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

cu 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 


PROGRAMMER AND DATE 

6. L. HOUSER 
Tkw SYSTEMS 
OUNt 1972 


PDKPOSt 

Subroutine, ploter is used to draw straight lines from a 

BASE POSITION TO EACH SCALED Y VALUE IN THE INPUT DATA 
ARRAY. 

USAGE 

Call PLOTER CNPTS) 

DESCRIPTION of parameters 
Input 

CALLING SEQUENCE 

NpTS - NUMBER OF data POINTS C STRAIGHT LINES) TO BE 
PLOTTED 

COMMON 

COMMON /A/ DAT AC 480) - ARRAY OF Y DATA POINTS 

TO BE PLOTTED. 

OUTPUT 

NONE 

REMARKS AND RESTRICTIONS 

DATA POINTS IN DATA ARRAY MUST HAVE VALUE BETWEEN 
1. AND 8. NO ERROR MESSAGE IS GIVEN. 

SUBROUTINES REOUIRtD 
Pt-N 

PLOTPT 


METhOD 

EACH DATA VALUE IS COMPARED TO THE BASE POSITION* IF IT 
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lb E«UmL» ThE X POSITION COUNTER IS INCREMENTED# IF THt. 

value is oreater* a line is drawn from the base position 

T® fHc Y uAT A VALUE# AND THEN THE POSITION COUNTER IS 
INCREMENTED* 

COMMON/ a/ tiATAmao) 

INTEGER UP* down 

npt=npts 

11=1 

DOwN=b 

uP=7 


DATA POINT VALUtS MUST BE BETWEEN 1» AND 8. 


YB1AS=1.0 

ASTART=**0.*i>2 

UX=.G333 


SELECT PEn 1# PLACE PEN IN UP POSITION 


CALL PENCttl 
CALL PENiUPJ 


MOVE PEN TO START POSITION 


CALL PLCTPH*START»YBIAS> 


EXAMINE AND PLOT EACH INPUT POINT 


X=XSTART 
DO SO I=U«FT 
t=uATA£W 

lF£Y-YBl«9fr ao*30»20 
20 CONTINUE 

CALL PlOTPHX* YblAS) 
CALL PENCDOW4) 

CALL PLOTPlCX»Y) 

CALL PENCUP) 

30 CONTINUE 




X=X+UX 
bU CONTINUE 


PLTGEN PROGRAM LISTING 


L 

c place pen in up position and return 

c 

CALL PEMUP) 

KLIURN 

tND 


3-35 



PLTGEN PROGRAM LISTING 


SUbKOUT I NE NUMdtRl XNUM t IFORM t X X # Y Y ) 

CL* ********************************** 


CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CO 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CO 

CD 

CD 

CD 

CD 

CD 

CD 

CD 


PROGRAMMER and date 

6* L» HOUSER 
TR# SYSTEMS 
JDNE 1972 


PURPOSE 

SUbROUTlNE NUMbER CONVERTS A FLOATING POINT NUMBER INTO 
A FORMATTED HOLLERITH STRING. AVAILABLE SPECIFIED FORMATS 
ARE F7.1 AND E1Q.3. 

USAGE 

CALL NUMBER t XUUM* IFORM* XX» YY ) 

DESCkIPTION OF PARAMETERS 


iNPDl 

CALLING SEUUENCE 

XNUM - FLOATING POINT NUMBER TO BE CONVERTED 

IFORM - FORMAT SELECTOR 0 = F7.i» 1 2 £10.3 

X - X COORDINATE OF FIRST DIGIT 

Y - Y COORDINATE OF FIRST DIGIT 


OUTPUT 

NONE 


REMARKS AND RESTRICTIONS 

SUBROUTINE NUMBER DOES NOT CHECK FOR* OR PROTECT AGAINST 
THE INPUT NUMBER EXCEEDING THE SPECIFIED FORMAT LIMITS. 

SUBROUTINES REQUIRED 
PRINT 


MEIHuD 

the sign and exponential magnitude of the input number 

IS DETERMINED. EACH DIGIT OF THE NUMBER IS THEN ISOLATED 
AND CONVERTED TO HOLLERITH DATA. THE PROPERLY FORMATTED 
NUMBER CONVERTED TO A HOLLERITH STRING IS THEN SENT TO 
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CD SUbKOuTINL PRINT. 

CL) 

CD* ********************************** 
DIMENSION ICODEt 14) » ICHARt 10 ) » IFIELDl 80 ) 

DATA LICODLLI) »I=lfl4)/5H0 »5Hl »5H2 »5H3 ,5H4 

* t SH5 »brtb * SH7 rSH8 »5H9 »5H+ »5H- »5H. 

* » ShE / 

DATA XbLANk/5H / 

10=0 
11=1 
12=2 
17=7 
110=10 
1=0 

NPOIm.R=u 
lPO*ER=0 
DEeX=.0666 
x=ax 
Y=YY 
10k 1N=1 
1FRMAT=1F0KM 
XN=XNUM 
DO 5 K=1 » 80 
IF IELDC K ) = it)LANK 
b CONTINUE 

DETERMINE IF NUMBER IS NEGATIVE AND SET SIGN WORD 

1S=IBLANK 

1FLXN-0.0) 30 » 100* SO 
30 CONTINUE 

1S=1 CODEC. 12) 

SO CONTINUE 

xn=abslxn) 

DETERMINE IF NUMBER HAS A NEGATIVE EXPONENT AND SET 
EXPONENT SIGN WORD ... 

IFLXN-1.0) 60'80'8Q 
60 CONTINUE 
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ISP-J CvuCt i<d) 

CALCULATE EXPONENT IF IT IS NEGATIVE 


6b 1=1+1 

APO*LK=iO.**I 
POfttK=l«/XPO«CK 
IFL AN+l^E-^O-PUrtEK) 65»70»70 
7U CONTINUE 
IPOhLR=i 
NPu«£R=iP0«ER+4 
GO TO 100 

bcT POSITIVE EXPONENT SIGN WORD AND DETERMINE EXPONENT 
VALUE 

O0 CONTINUE 

iSP=lCou£C 11 ) 
e>5 1=1+1 

POwEK=10.**I 

iFLXN+l.E-o-POwER) 90 r 85 t 65 
90 CON1INUE 
lPO*tK=l“l 
UPOwEK=*t-IPOrtER 
100 CONTINUE 

iX=XN*10.**NPOrtLR+l.E-3 

ISOLATE DIGITS OF NUMBER BEING CONVERTED 

uO 110 0=1 » 5 
K=0” J 

1TEMP=IX/10 

N=IX-ITEMP*10+1 

1CHARCK)=IC0DELN) 

IX=1TEMP 
110 CONTINUE 


IF F7.1 FORMAT GO TO 120 
1FLIFRMAT.EO.0) GO TO 120 
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SET UP HOLLERITH STRING FOR E1Q.3 FORMATTED NUMBER 

1TcMP=1P0WER/1U 

N=IPOw£R-ITEMP*iO+l 

1CHARE6)=1COQECn) 

1X=ITEmP 
1TEMP=IX/10 
h= 1X-ITEMP*10+1 
ICHAKCb) =1 CODELh) 

1FIELDL1)=IS 

aFIELD12)=1CHAKL1) 

IF ICLUL 3) =1CUUCL 13) 

IF ltLUL i +) = lCHARC.2) 

1F1EluL5> = 1CHMKC3) 

IFIEluL 6 ) =IChMRL 4 ) 

IF IElUL 7 ) =ICOUEL 1H ) 

lFlELblo)=iSP 

IF IELUL 9 ) =alhAKL 5 ) 

1FIElUL10)=ICHAKCo) 

CALL PrtiNTt X# Y# IFIELbt I10» I OR IN »DELX ) 

RETURN 


SET UP HOLLERITH STRING FOR F7.1 FORMATTED NUMBER 


12U CONTINUE 

1ST=4-1P0*CR 

u=0 

uO 160 1=1 » 7 
IFC I~IST ) 16Q» 130 » 1A0 
130 CONTINUE 

IFlELDt 1)=1S 
GO TO 160 
1HU CONTINUE 

IFCl.NE.fa) GO TO 150 
IFIElDC b) = ICOOEt 13) 
GO TO lbO 
ISO CONTINUE 
J=J+1 

IFIELUC1)=ICHARCJ) 
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ibb CONTINUE. 

t ACL Pr\ i NT t X » Y » I F I ELD » 1 7 . 1 OR I N » DELX > 

RETURN 

END 


V 



uu« * 

CU 

Cu 

CU 

cu 

cu 

cu 

Cu 

Cu 

CO 

Cu 

Cu 

cu 

CU 

CO 

CO 

Cu 

cu 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CD 

CO 

CO 
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SUBROUTINE PL01PT IX»Y) 

********************************* 

PROGRAMMER AND DATE 

G. L. HOUSER 
TRW SYSTEMS 
June i.y72 

PURPOSE 

subroutine plotpt is used to give the pen a command to 

MOVE FROM ITS PRESENT POSITION TO A NEW POSITION CX»Y). 
The PEN MAY BE IN EITHER THE UP OR DOWN POSITION. 

oSmGE 

CmLL PlOThT CXrY) 
uESCRlPTIOu uF PARAMETERS 
INPU1 

CALLING sequence 

X - mOSISSA VALUE BETWEEN -9,999 AND +9.999 

Y - OROINATE VALUE BETWEEN -9.999 AND +9.999 

OUTPUT 

NONE 

kEmARKS ANu RESTRICTIONS 

IF THE X OR Y VALUES EXCEED THE ABOVE LIMITATIONS THE 
PLOTTER WILL RETURN TO THE CENTER OF THE PLOT BOARD AND 
CONTINUE PLOTTING. 

SUBROUTINES REQUIRED 
COMAND 

METHOD 

THE X AND Y COORDINATES ARE SCALED AND SENT TO SUBROUTINE 
COMAND. 



CD* ********************************** 
11=1 





PLTGLN PkOGkAM LISTING 


XX=X 

YY-Y 

lX=XX*lUGO. 

iY=YY*lUOU. 

CALL CGKANuL il * i.1 r 1X» XY ) 

KLlUk.x 

LNL 



c 


o 


PLTGEN PROGRAM LISTING 

SUuKUOUNE InITAL L IaSCAL’ IYSCAL# IXLOS’ IYDOS* IXBOS* IYBOS) 


CL 



CL 

PROGRAMMER AWL UATE 

cu 


G. L. HOUSER 

CL 


Tk* systems 

CL 


June 1972 

CL 



CL 

PURPOSl 


CL 


THE PuRPOSc OF SUBROUTINE INITAL IS TO INIT T ALlZE 

CL 


THE PLOT PROGRAM BY SETTING UP SCALE FACTORS* "’OARD 

CL 


OFFSET » Anu LATA OFFSET COMMANDS. 

CL 



CL 

LSAGc 


Cu 


CALL InITAL £IXSCAL*IYSCAl»IXDOS* IYDOS* IXBOS* IYBOS) 

CL 



CL 

lEs>Cr!PT ION OF PARhMEIERS 

CL 



CL 

Input 


CL 

CALLING SEGUEnCE 

CL 


IaSCAl - X SCALE FACTOR £3333 IS NOMINAL ONE TO ONE 

Cl 


CORRESPONDENCE) 

Cl 


1YSCAL - Y SCALE FACTOR tSAME AS ABOVE) 

CL 


1XDOS - X LATA OFFSET £0000 IS NO OFFSET) 

CL 


lYuOS - Y LATA OFFSET C0O00 IS NO OFFSET) 

Cu 


IXBOS - X BOARD OFFSET CSAME AS ABOVE) 

CL 


IYdOS - Y BOARD OFFSET CSAME AS ABOVE) 

CL 



CL 



CO 

OUTPUT 

CL 


NONE 

CL 



CL 

REMARKS 

AND RESTRICTIONS 

CL 


NONE 

CL 



CL 

SUBROUTINES RtaUIRLL 

Cl 


COMANO 

CL 



CL 

METHOD 


CU 


THE SCALE FACTORS# BOARO OFFSET #AnD DATA OFFSET 


* * * 


£ 


C 
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CD commands are formatted and sent to subroutine .comano. 

CD 

CD* *%*****♦*******<•************♦***** 

DIMENSION NW0KuL2) 

IXS=l asCAL 

XYS-IYSCAL 

IXO=1XDOS 

IYD=1Yuw»S 

IXD=IXbuS 

lYtPIYouS 

il=l 

12=2 

14=4 


SET UP SCALE FACTOR 


aSC=2 

CAtL COMANDL ISC» ISC» IXS» IYS) 
SET UP DATA OFFSET 


1D0C=3 

CALL CohANOC IDOC • IDOC » IXU» I YD) 
SET UP BOARD OFFSET 


1B0C=4 

CALL COMANdC IDOC t IbOC * IXB» I YB) 

RETURN 

END 
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CU* 

CO 

CD 

cu 

Cu 

Cu 

cu 

CD 

CU 

cu 

CU 

uU 

CU 

uu 

Cu 

Cu 

Cu 

CU 

CU 

cu 

CU 

cu 

Cu 

Cu 

Cu 

Cu 

Cu 

cu 

cu 

Cu 

Cu 

CU 

cu 

cu 

Cu 

CU 

cu 

Cu 

cu 


SUUKUL I i HE PLNt 1 SLLC T ) 

********************************** 

PKOURAlwtK AnU UATl 

G. L. HOUSER 
lit# SYSTEwb 
JUNE IV72 


PURPOSL 


USAUL 


SUUROcIT 1NE PEN IS USED TO SELECT EITHER THE PLOT PEN OR 
THE PR INTER » AND TO GIVE PEN UP OR PEN DOWN COMMANDS. 


CAuL PtN LISLECT) 


uESCkIPT 1<Ji4 ul- PARAMETERS 

Input 

calling sequence 

ISLECT - 1 PEN 1 C PLOTTER) 
2 PEN 2 C PRINTER) 
b PEN Down 
7 PEN UP 

OUTPUT 

NONE 

REMARKS ANu RESTRICTIONS 

none 

aUoROU TINES Rt GUI RED 
COMAND 


rtEThOD 


ISLECT IS CHECKED TO DETERMINE IF HE REQUEST IS A SELECT 
PEN OR PRINTER COMMANQ# THE APPROPRIATE COOE IS SENT 
TO SUBROUTINE COMAND# ISLECT : o 0? 7 AUTOMATICALLY C-IVES 
PEN UP OR PEN DOWN COOE 




& 

I 


r 
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Cl>* ********************************* 
iS=ISLtCT 
1X-U 
iY=ls 

IS = ir SELECT PEN 1 WITH CODE 11 
lFllS.Ev».l> IY=11 

IS = 2f SELECT PRINTER wITH CODE 12 

lFLU.Eu.2) iY-12 
C 

CALL COhANuClX'iWIXtIX) 

C 

KtTUKN 

tNU 
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SUoKOU l i l.E PR I N 1 L X * T * ISTRNG » NUMbER * I OR I N * DEL ) 


Lb* 

***** 

* + *************************** 

CU 



Lb 

PhOGhAMi v *ER ANb bATL 

Cb 


G. L. HOUStR 

Lb 


Trw systems 

Cb 


UunE 1972 

Cb 



Cb 

PUKPOBl 


CO 


SUBROUTINE PRINT IS USED TO CONVERT A STRING OF HOLLERITH 

Cb 


characters into the eat. plotter character set. the 

Cb 


CONVERTED CHARACTERS ARE then SENT TO SUBROUTINE COMANO 

Cb 


To bE FORi v iEb InTO PLOT COMMANDS. 

CU 



Cb 

USAGE 


Cb 


CALL PRINT L X* Y * ISTRNG » NUMBER » IORIN»DEL) 

Cb 



cu 

INSCRIPTION OF PARmMETERS 

CU 



cu 

I.NiPUt 


Cb 

CALLING SEQUENCE 

CU 


X - X COORDINATE VALUE* STARTING POSITION BETWEEN 

Cb 


-9.999 AND +9.999 

Cb 


Y - Y COORDINATE VALUE* STARTING POSITION BETWEEN 

CU 


-9.999 AND +9.999 

Cb 


ISTRNG - ARRAY CONTAINING CHARACTER STRING 

cu- 


NUMBER - NUMBER OF CHARACTERS IN CHARACTER STRING 

cu 


IORIN - ORINTATION OF CHARACTERS 0 - VERTICAL* 

cu 


1 - HORIZONTAL 

cu 


DEL - DISTANCE BETWEEN CHARACTERS IN INCHES 

cu 



Cb 

OUTPUT 

Cb 


NONE 

Cb 


- 

CU 

REMARK^ 

ANU RESTRICTIONS 

Cb 


The MAXIMUM NUMBER OF CHARACTERS TO BE CONVERTED IN ONE 

cu 


STRING IS 80. 

cb 



Cb 


The CHARACTER STRING MUST BE CONTAINED IN AN ARRAY HAVING 

CD 


ONE CHARACTER PER WORD 
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w l»Ut>ROUTiNES REOUIREu 

.J Com AND 

. j 

v. ML Thou 

each character in the string is compared with each 

•j CHARACTER In THE EAI CHARACTER SET. WHEN a MATCH IS MmJE 

Cb THE CODE IS SENT TO SUBROUTINE COmAND. THE SPACING 

-t- VARIABLE IS INCREMENTED AND THE NEXT CHARACTER IN THE 

SlRINv? IS PROCESSED. 

uL> 

.* ********************************** 

DIMENSION ISTRNGCBU) . ICHARC4S) » IC0DEC4S) .LABELC80) 


EhI character table 


UATA 

t 

ICHAKC I 

) »I=1»44)/5H0 

• 5h1 

» 5h2 

1 5H3 

»5H4 

9 

♦ 


bHb 

» brto 

rbH7 

»5H8 

» 5H9 

1 5HA 

1 5HB 

9 

* 


bhC 

#SHD 

.SHE 

»5hF 

1 5HG 

»5HH 

.5Hl 

9 

* 


5HJ 

» bHK 

»bHL 

»5HM 

1 5HN 

1 5H0 

»5HP 

9 

* 


ShU 

» bHR 

» bHS 

»5HT 

1 5HU 

»5HV 

»5HW 

9 

* 


5HX 

» t>HY 

»5HZ 

*5H- 

»5H+ 

»5HC 

1 5H) 

9 

♦ 


bH= 

» bH* 

»5H» 

»5H. 

/ 




UATA 

C 

1CODECI 

).1=1»44)/ 46. 

it 2. 3 » 

4. 5. 

6 » 7. 10. 11 .61 .62 » 

63 

* .64 

»b5» 66 » 

b7» 70 » 71 

.41 »42i 

>43»44»45i 

.46.47. 

50»51»22» 

23.24.25 

fZb 


* * 27. 30 . 31 » 4u » bl) . 34. 74 »13» 54.33*73/ 

C 

DATA I6LANK/5H / 

NUM-NUMbER 
1D=10RIw 
DELT A=ucL 
10=0 
11=1 

lSLECT=b 

XX=X 

YY=Y 

uO 10 1=1 »NUM 
lABELC I ) =1STRNGC I ) 

10 CONTINUE 
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IDENTIFY EACH CHARACTER I.m THE STRING 

uO bo 1=1 »nUM 

1X=XX*100Q. 

iY=YY*luuO. 

IF CHARACTER IS A BLANK » INCREMENT CHARACTER SPACER 

IFCLAbELL 1 ) »£tt* 1 BLANK) GO TO 50 

SEARCH EAI CHARACTER TABLE 

DO 3u o=i » 44 
N=J 

lFi.LAbc.i~L 1 ) «E0» aCHAKCn) ) GO TO 40 
50 CCul xwuE 
40 CONTINUE 


load character and orientation code 
ich=icodecn) 

IFL lU.Eu.l) ICH=-ICH 

SEND SELECT CHARACTER COMMAND 
CALL COHANuC IQ » ISLECT » 10 » ICH) 

send plot character command 

CALL COMANOCll'll'iX'IY) 

60 CONTINUE 

iFLlu.EG.l) GO TO 55 

Increment y spacer if horizontal letters c vertical strng) 

YY=Y Y+utLTA 

GO TO bO .. • 

bo CONl INUE 


INCREMENT X SPACER IF VERTICAL LETTERS C HORIZONTAL STRNG) 
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AX— XX+L/t-t-T A 
oil lOi* I iiiJC. 
t 

i«LTU\N 

LNl) 



C 
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^UoKguUimL MAGTAPLo) 

S.U* ********************************** 

4 


W U 

KKObKMi v ti*iLK A NO 

LATl 



V-U 

KaChAKL C. THOMS 




TRW SYSTf.'iS 




Jo.C 1972 



CD 





.u 

PlkPosl. 




U-J 

MAGTAP 

IS M DUMMY ROUTINE that ALLOCATES CORE STORAGE FOR i 


V 

A CS-1 

ASSEMBLY LANGUAGE ROUTINE NAMED MAGTAP. 

THE CS-1 


' 

PkOGM 

is PuFORMS ALL TAPE OPERATIONS REQUIRED BY 

the pro- 


CL 

gram* 

hUK Tape parity ERRORS the MAGTAP ROUTINE 

attempts 



KLCOVEkY M maximum of ten times 




FwK UNktCOVtKABLE PARITY uRKORS* FRAME COUNT ERRORS ANu 1 



T XkINO 

ERRORS THE COMPUTER COMES TO A HALT WITH 

NO PRINT- 


LG 

Eu EKRwK .-.ESSAGES. THE OPERATOR IS ABLE TO DETERMINE THE 1 


LU 

NATURE 

OF THE ERROR BY INDICATOR LIGHTS ON THE 

TAPE DRIVE 


CL* 





CL 

USAgc. 




CU 

CauL MAGTAPCIOPT»NUM» ARRAY* IUNT) 



Cu 





CL 

uESCrtlPTION OF 

PARAMETERS 



CU 





CU 

INPul . 



<*' 

CD 

CALL 1 bduUENCE. 



Cu 

lopT 

- »EQ » 1 READ TAPE RECORD 



LG 

IOPT 

- »t.Q»2 WRITE TAPE RECORD 



LU 

I OPT 

- .EQ.3 BACKSPACE RECORD 


*)S. - 

CU 

IOPT 

- .tQ.4 REvtlND TAPE 



CD 

lop r 

- .EQ.5 WRITE FILE MARK 


fa 

CD 

IOPT 

- .tQ.6 SPACE FILE FORWARD 



CU 

IOPT 

- .EQ.7 SPACE FILE BACKWARD 



CD 

NUM 

- NUMBER OF WORDS IN DATA RECORD 


r 

CD 

ARRAY 

- DATA RECORD 



CD 

IUNT 

- TAPE UNIT 



CD 

TAPE 


! 


CU 





CD 

OUTPUT 




CD 

Tape 


♦ 

- 
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Cu 
CD 
C 0 
Ou 
CL) 
CD 
CO 
Cu 
cu 
CO 
cu* 


KLmARNS AND RESTRICTIONS 

WhcN TaPE OPERATION IS NEITHER A READ OR WRITE ». NUM ANj 
AkraY ARE DUMMY parameters. 

This routine allocates the necessary core required for 

THE CS-1 PROGRAM aHICH MUST bE LOADED BY PAPER TAPE. SEE 
PauoKAh OPERATING INSTRUCTIONS FOR CORRECT PROCEDURES. 
SDbKwU TINES REQUIRED 
none 

********************************* 

DIMENSION mL300) 
aLD = 1.0 
re IUKn 
Enu 
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